(19) 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(id EP 1 043 687 B1 

EUROPEAN PATENT SPECIFICATION 



(45) Date of publication and mention 
of the grant of the patent: 
22.11.2006 Bulletin 2006/47 

(21 ) Application number: 00400949.4 

(22) Date of filing: 06.04.2000 



(51) IntCI.: 

G06T 1/00(2006.01) 



(54) Method for inserting a watermark and associated decoding method 

Verfahren zum Einfugen eines Wasserzeichens und dazugehoriges Dekodierverfahren 
Methode pour inserer un filigranne et methode associee de decodage 



CD 

CO 
CD 

CO 



CL 
LU 



(84) Designated Contracting States: 
DE FR GB 

(30) Priority: 09.04.1999 FR 9904461 
09.04.1999 FR 9904462 
25.06.1999 FR 9908187 

(43) Date of publication of application: 
1 1 .1 0.2000 Bulletin 2000/41 

(73) Proprietor: CANON KABUSHIKI KAISHA 
Tokyo (JP) 

(72) Inventors: 

• Donescu, loana 
35000 Rennes (FR) 

• Nguyen, Eric 
35000 Rennes (FR) 

• Onno, Patrice 
35000 Rennes (FR) 

(74) Representative: Santarelli 

14, avenue de la Grande Armee, 
B.P. 237 

75822 Paris Cedex 17 (FR) 



(56) References cited: 

• MENG AND CHANG: "EMBEDDING VISIBLE 
VIDEO WATERMARKS IN THE COMPRESSED 
DOMAIN" INTERNATIONAL CONFERENCE ON 
IMAGE PROCESSING ICIP98, vol. 1,4-7 October 
1998, pages 474-477, XP0021 24674 NEW-YORK 
NY US 

• SWANSON E.A.: "OBJECT-BASED 
TRANSPARENT VIDEO WATERMARKING" IEEE 
FIRST WORKSHOP ON MULTIMEDIA SIGNAL 
PROCESSING, 23 - 25 June 1997, pages 369-374, 
XP0021 24675 NEW-YORK NY US 

• PODILCHUK C I ET AL: "IMAGE-ADAPTIVE 
WATERMARKING USING VISUAL MODELS" IEEE 
JOURNAL ON SELECTED AREAS IN 
COMMUNICATIONS,US,IEEE INC. NEW YORK, 
vol. 16, no. 4, 1 May 1998 (1998-05-01), pages 
525-539, XP000765113 ISSN: 0733-8716 

• TAO B ET AL: "ADAPTIVE WATERMARKING IN 
THE DCT DOMAIN" IEEE INTERNATIONAL 
CONFERENCE ON ACOUSTICS, SPEECH, AND 
SIGNAL PROCESSING (ICASSP),US,LOS 
ALAMITOS,CA: IEEE COMP. SOC. PRESS, 1997, 
pages 2985-2988, XP000788033 ISBN: 
0-8186-7920-4 

• RUANAIDH J J K O ET AL: "WATERMARKING 
DIGITAL IMAGES FOR COPYRIGHT 
PROTECTION" IEE PROCEEDINGS: VISION, 
IMAGE AND SIGNAL PROCESSING,GB, 
INSTITUTION OF ELECTRICAL ENGINEERS, vol. 
143, no. 4, 1 August 1996 (1996-08-01), pages 
250-256, XP000613938 ISSN: 1350-245X 



Note: Within nine months from the publication of the mention of the grant of the European patent, any person may give 
notice to the European Patent Office of opposition to the European patent granted. Notice of opposition shall be filed in 
a written reasoned statement. It shall not be deemed to have been filed until the opposition fee has been paid. (Art. 
99(1) European Patent Convention). 



Printed by Jouve, 75001 PARIS (FR) 



EP 1 043 687 B1 



Description 

[0001] The present invention concerns a method of inserting a watermark and an associated method of decoding a 
watermarking signal. 

5 [0002] Correlatively, it concerns a device for inserting a watermarking signal and a device for decoding a watermarking 
signal. 

[0003] The present invention lies in general terms in the technical field of watermarking of digital images, and more 
particularly still images. 

[0004] Watermarking digital data makes it possible to protect these data, for example by associating copyright infor- 
10 mation with them. 

[0005] In its general principe, watermarking consists of inserting an indelible watermark in digital data, which can be 
likened to encoding additional information in the data. 

[0006] Decoding this additional information makes it possible to check the additional information which has been 
inserted. 

15 [0007] This inserted watermark must consequently be both imperceptible, robust to certain distortions applied to the 
digital image and reliable to detect. 

[0008] Conventionally, a usual technique of inserting a watermarking signal in a digital image consists of using a linear 
modulation model in which at least one subset of coefficients representing the digital image is modulated according to 
this linear model using a weighting coefficient. 
20 [0009] By denoting as X = {Xj, 1 < i < N} a set of the coefficients representative at least of part of a digital image, and 
as w = {Wj, 1 < j < P} a watermark of size P < N, a pseudo-random signal with known distribution and average of zero, 
the linear insertion formula is : 



25 X^Xj+bajWj.with 1 <j<P, 

in which {Xj, 1 < j < P} i s a subset of the set of coefficients X, b is an information bit, and ocj is a weighting coefficient, 

also called modulation amplitude. 
30 [001 0] The detection of the watermark then consists in detecting whether or not the pseudo-random sequence w has 

been inserted into a set of coefficients. This detection is carried out without using the original image and can be based 

on a standardised statistical test which makes it possible to calculate a probability of correct detection. 

[0011] Such an insertion technique makes it possible, by insertion of a watermark, to insert a single information bit 

since the response from the detector is binary (yes/no). 
35 [0012] In order to insert a larger number of information bits into the digital image, in particular when a code of Q bits 

is wanted, indicating, for example, the name or the address of the owner or of the author of the image, it is necessary 

to reiterate the previously described insertion method as many times as there are information bits to be inserted. Typically, 

in order to insert a binary signal, either b = 1 , or b = -1 is used. 

[0013] Put another way, Q subsets of coefficients have to be chosen and the modulation of these subsets has to be 
40 carried out by choosing Q watermarks. 

[0014] Separate subsets of coefficients are preferably chosen such that the modulation operations are not superim- 
posed on one another, which could disturb the detection or cause troublesome visual effects. 

[0015] It is a matter, consequently, of choosing a partition of the coefficients representative of the digital image into 
Q separate subsets, each carrying one information bit. 

45 [0016] Numerous known methods use a technique of inserting a watermarking signal of given size by spectrum 
spreading. The disadvantage of these methods, describe for example in the article entitled "Secure spread spectrum 
watermarking for multimedia" by I.J. COX et al, in Proc. ICIP, pages 243-246, Sept. 1996 and in the article entitled 
"Digital watermarking of raw and compressed video" by F. HARTUNG et al, in Proc. SPIE 2952 : Digital Compression 
Technologies and Systems for Video Communication, pages 205-213, October 1996, is that they use an arbitrary 

50 partitioning of the image into blocks of fixed size without any guarantee of detectability of the modulations effected on 
each of the blocks. 

[0017] The aim of the present invention is to propose a method of inserting a watermark wherein a partitioning of a 
digital image is adapted to the image, ensuring detectability of the bits inserted for a technique of insertion by modulation. 
[001 8] ' The present invention thus relates to a method of inserting awatermarking signal in adigital image by modulation 
55 of coefficients of distinct regions representing said digital image, comprising a step of determining a partitioning into 
distinct regions of the digital image. 

[0019] In accordance with the invention, this partitioning into regions is effected by an adaptive partitioning as a function 
of a criterion of statistical detectability of an information bit inserted in each region. 
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[0020] In this way an adaptive partitioning is performed which has the advantage of adapting to the content of the 
image signal and guaranteeing detectability of the inserted watermarking signal, unlike the choice of an arbitrary parti- 
tioning into fixed blocks. 

[0021] This spatial adaptivity is necessary so as the content of the image is not spatially homogenous. 
5 [0022] This is because, for a given modulation amplitude, the greater the length of the sequence of the watermarking 
signal inserted, and consequently the greater the size of the region and therefore the number of modulatable coefficients, 
the more reliable is the detection. 

[0023] It can in fact be shown that, when a correlation calculation is used for effecting the detection, the response of 
the detector is proportional to P 1/2 where P is the length of the sequence of the watermarking signal. 
10 [0024] Moreover, the response of the detector is also dependent on the image signal itself. The greater the variance 
in the coefficients representing the image, the weaker the response of the detector. 

[0025] In conclusion, for a given modulation amplitude and image, there exists a minimum length of the watermarking 
signal necessary for ensuring a given detection probability, this minimum length also depending on the local content of 
the image signal in the region in which the watermarking signal is inserted. 
15 [0026] The adaptive partitioning effected by the method according to the invention guarantees the detectability of the 
watermarking signal inserted in the different regions determined by this partitioning. 

[0027] Preferably, the detectability criterion is calculated using a maximum value of the weighting coefficient ensuring 
imperceptibility of the watermarking signal for the modulation of the set of coefficients representing said digital image. 
[0028] This is because, the greater the modulation amplitude, that is to say the value of the weighting coefficient, the 
20 more reliable is the detection in that a better probability of detection can be obtained. Nevertheless, the greater this 
modulation amplitude, the more this modulation may be visible in the watermarked image. 

[0029] It is therefore judicious to choose a modulation amplitude equal to a maximum amplitude, generally denoted 

JND (Just Noticeable Difference), beyond which the modulation becomes visibly perceptible in the image. 

[0030] It can be shown that this maximum amplitude also depends on the length of the sequence of the watermarking 

25 signal and that it decreases with this length. 

[0031 ] By choosing a maximum value of the weighting coefficient ensuring imperceptibility of the watermarking signal 
for the modulation of the set of coefficients representing said digital image, the value of the weighting coefficient on each 
region of the partitioning, and with a size less than the cardinal number of all the coefficients representing the image, 
will be much less than the maximum modulation amplitude which guarantees imperceptibility of the modulation on this 

30 region. 

[0032] The partitioning thus performed also ensures invisibility of the modulation effected on each of the regions 
obtained by the partitioning. 

[0033] According to an advantageous characteristic of the invention, the detectability criterion is calculated using, for 
each coefficient to be modulated, a weighting coefficient according to a law of the form <Xj = kj.oty, where kj is a modulation 
35 factor dependent on the coefficients situated in a vicinity of the coefficient considered on said region and oc v is equal to 
said maximum value of the weighting coefficient. 

[0034] It is thus possible to exploit the possibilities of visual "masking" of the modulation by the image signal itself. 
Typically, use is made of the fact that, the greater the activity of the image signal, in the sense for example of a significant 
local variance in the coefficients, the less visible will be the modulation. 
40 [0035] It is thus possible to locally increase the weighting coefficient and therefore the modulation amplitude, to the 
benefit of detectability, without impairing the imperceptibility of the inserted modulation. 

[0036] According to another preferred characteristic of the invention, the method of inserting a watermark includes a 
step of applying a distortion before calculating the detectability criterion. 

[0037] According to another preferred characteristic of the invention, the method of inserting a watermark includes a 
^5 prior step of comparing the size of each region of the partitioning with a minimum size corresponding to the minimum 

size of a statistically significant sample for detecting a watermarking signal inserted in said region. 

[0038] Indeed, because of the calculations of a statistical nature performed when the watermarking signal is detected, 

notably when standardised hypothesis tests are used, it is necessary to consider a minimum boundary for the size of 

each region of the partitioning guaranteeing the validity of the detection tests. 
50 [0039] According to one particularly practical embodiment of the invention, for each region of the partitioning, all the 

coefficients of said region are modulated by inserting the watermarking signal, a detectability amplitude is calculated 

from said detectability criterion and said detectability criterion is validated by comparing said detectability amplitude with 

a predetermined threshold value. 

[0040] The detectability amplitude can be the result of a calculation of correlation between the set of modulated 
55 coefficients and the watermarking signal or the result of a standardised statistical test such as the one described in the 
article entitled "A method for signature casting on digital images" by I. PITAS, in Proc. ICIP, pages 215-218, September 
1996. 

[0041] Practically, the detectability criterion can be validated if a minimum modulation length calculated for ensuring 



3 



EP 1 043 687 B1 



detectability of the inserted watermarking signal in less than or equal to the cardinal number of said region. 
[0042] According to a preferred embodiment of the invention, the method comprises the following steps: 

generation of a centred pseudo-random sequence equal in size to the cardinal number of a set of coefficients 
5 representative of at least a part of the image, formed from centred pseudo-random sub-sequences; 

modulation of the said set of coefficients by the said centred pseudo-random sequence in order to insert the same 
information bit on the said set of coefficients; and 

checking, over each region, of a criterion of detectability of the information bit inserted by modulation. 

10 [0043] It is thus possible to achieve a partition of a set of coefficients into distinct regions, by carrying out an adaptive 
partitioning on the basis of a criterion of detectability of an information bit inserted on each region. 
[0044] Such an adaptive partitioning exhibits the advantage of adapting to the content of the image signal and of 
guaranteeing the detectability of the watermark inserted on the different regions determined by this partition, in contrast 
to the choice of a random partitioning. 

15 [0045] In practice, for each region of the partition, a criterion of detectability of an information bit inserted by modulation 
of the coefficients of the said region is checked. This criterion may consist in the comparison of a result of a standardised 
statistical test with a threshold value corresponding to a probability of detection fixed in advance. 
[0046] By virtue of the method in accordance with the invention, it is possible to generate a centred pseudo-random 
sequence on the basis of the merging of centred pseudo-random sequences of smaller sizes, and to modulate all the 

20 available coefficients in order to achieve the insertion of a watermark. 

[0047] By modulating all the coefficients with the same information bit and a centred pseudo-random sequence which 
is broken down into pseudo-random sub-sequences which are also centred on each region of the partition produced 
over the set of coefficients, the steps of generating a pseudo-random sequence and of modulation do not have to be 
reiterated at each new region tested. 

25 [0048] According to one preferred characteristic of the invention, the method of inserting a watermark comprises a 
single step of applying a distortion to the said set of modulated coefficients, before the step of checking of a detectability 
criterion. 

[0049] It is known to apply a predetermined distortion to a set of modulated coefficients and to carry out the check on 
the detectability criterion after this pre-distortion, in such a way that it is possible to guarantee the detection of the 
30 watermark inserted on a partition of the coefficients for a predetermined distortion. 

[0050] By virtue of the invention, a pre-distortion can be applied to the set of modulated coefficients, in a unique way, 
and does not need to be reiterated at each new region tested in the partition. 

[0051] In one particularly practical embodiment of the invention, the set of coefficients is a set of spectral coefficients 
originating from a spectral transformation of a digital image, and the step of applying a distortion includes the following 
35 sub-steps: 

inverse spectral transformation of the said set of modulated coefficients in order to obtain a watermarked image; 
application of a distortion to the said watermarked image; and 

spectral transformation of the watermarked image in order to obtain a set of modulated coefficients after distortion. 

40 

[0052] The saving obtained in terms of complexity of calculation is all greater in this case whenever the application of 
a distortion requires several steps of processing of the set of coefficients, and particularly the application of an inverse 
spectral transformation in order to obtain the watermarked image and to apply a distortion, such as compression/de- 
compression to it, and then the application of a forward spectral transformation of the image again so as to check the 
45 criterion of detectability in the spectral domain. 

[0053] According to another particularly practical characteristic of the invention, the centred pseudo-random sequence 
is formed from pseudo-random sub-sequences of identical size. 

[0054] Such a pseudo-random sequence is particularly well adapted to the application of a bottom-up or top-down 
iterative partitioning of the quadtree type which makes it possible iteratively to combine regions of identical sizes or 
50 iteratively to divide regions into sub-regions of identical sizes. 

[0055] For preference, the pseudo-random sub-sequences have a size greater than or equal to a minimum number 
of coefficients forming a region of the said partition of a size which is statistically significant for validating the said criterion 
of detectability over the said region. 

[0056] It is thus possible theoretically to form a maximum partition of the set of coefficients with the same unique 
55 centred pseudo-random sequence, each pseudo-random sub-sequence being centred on the regions with the smallest 
size which is statistically admissible for inserting a watermark by modulation into the set of coefficients. 
[0057] According to one practical embodiment, which is inexpensive in terms of calculations, at the generation step, 
each centred pseudo-random sub-sequence is created by generating half of the pseudo-random samples of the said 
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sub-sequence via a pseudo-random number generator and by symmetrising the said pseudo-random samples generated. 
[0058] For preference, in this embodiment, at the generation step, each centred pseudo-random sub-sequence is 
created by moreover permutating the said pseudo-random samples using a secret key. 

[0059] This permutation operation makes it possible better to preserve the intrinsic random-sequence nature of a sub- 
5 sequence created by symmetrising pseudo-random samples, and not to diminish the imperceptibility of the inserted 
watermark. 

[0060] According to another practical embodiment of the invention, at the generation step, each centred pseudo- 
random sub-sequence is created by generating pseudo-random random samples via a pseudo-random number generator 
and by redistributing over each sample the sum of the said pseudo-random samples generated. 

10 [0061] According to a preferred characteristic of the invention, a partitioning by iterative division of the digital image 
is performed and, for each region at a rank of said partitioning, a partitioning of a directly higher rank of said region is 
validated if and only if there exist at least two subregions of said region for which said detectability criterion is validated. 
[0062] Thus, when it is attempted to divide a region into subregions, this partitioning of a higher rank is validated only 
if the number of regions in which a watermarking signal can actually be inserted whilst complying with the detectability 

15 criterion is actually increased. 

[0063] Otherwise, it is preferable to keep the region before division, of greater size, in order to increase the detectability 
of the watermarking signal. 

[0064] According to a preferred embodiment of the invention, for each region at a rank of said partitioning, the value 
of a capacity of the digital image is incremented when the partitioning of directly higher rank is not validated, the said 
20 capacity being equal to the cardinal number of the distinct regions obtained by the said partitioning on which the de- 
tectability criterion is validated. 

[0065] It is thus possible to calculate the capacity of the digital image which thus corresponds to the number of 
information bits which can be inserted in the digital image. 

[0066] By virtue of the adaptive partitioning performed by the invention, by reiterating the partitioning of each region 
25 as long as the detectability criterion is complied with, there is obtained the maximum capacity of the digital image which 
corresponds to the watermarking capacity of this image by the technique of insertion by modulation. This watermarking 
capacity can be defined as the maximum number of information bits which can be inserted imperceptibly in the image 
with a guarantee of detectability on decoding. 

[0067] According to another embodiment of the invention, for each region at a rank of said partitioning, the value of 
30 a capacity is incremented when the partitioning of directly higher rank is not validated, said value of the capacity is 
compared with a fixed capacity value and the partitioning is reiterated for a region having a detectability amplitude which 
is the h ighest amongst all the other regions to be processed when said capacity value is less than said fixed capacity value. 
[0068] This partitioning method makes it possible, in this embodiment, to find the best possible partitioning in terms 
of detectability when a predetermined number of information bits have to be inserted, equal to said fixed capacity value. 
35 [0069] By selecting, during the partitioning, the regions for which the detection amplitude is the highest, the present 
invention makes it possible to determine the partitioning which, a priori, guarantees the best detectability on decoding 
of the watermarking signals making it possible to insert a predetermined number of information bits. 
[0070] According to another embodiment of the invention, a partitioning is effected by iterative merging of the digital 
image and, for each region of said partitioning for which said detectability criterion is not validated, said region is merged 
40 with another region of said partitioning. 

[0071] Preferably, in this embodiment, said other region of the partitioning is if possible a region for which said de- 
tectability criterion is not validated. 

[0072] It is thus possible to increase the capacity of the digital image by joining the regions in which the detectability 
criterion is not satisfied. 

45 [0073] The present invention also concerns a method of decoding a watermarking signal in a digital image, inserted 
by modulation of coefficients representing said image in distinct regions forming a partitioning of said digital image. 
[0074] In accordance with the invention, this decoding method includes a step of determining a partitioning into distinct 
regions of the digital image to be decoded by effecting an adaptive partitioning as a function of a criterion of statistical 
detectability of an information bit obtained by demodulation of the coefficients on each region. 

50 [0075] Thus this decoding method makes it possible to find, at the decoder, the adaptive partitioning performed for 
inserting a watermarking signal. 

[0076] This is because, when the partitioning has been performed in an adaptive fashion as described above by a 
method of determining a partitioning in accordance with the first aspect of the invention, this partitioning, which depends 
on the content of the image, is not known to the decoder. 
55 [0077] The decoding method according to the invention then consists of testing all the acceptable partitionings and 
verifying the detection of a watermarking signal on each of the regions associated with this partitioning. 
[0078] According to a preferred characteristic of the invention, the decoding method includes a prior step of comparing 
the size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistically 
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significant sample for the detection of a watermarking signal inserted in said region. 

[0079] According to a practical embodiment of the decoding method according to the invention, for each region of the 
partitioning, a detection amplitude is calculated from said detectability criterion and said detectability criterion is validated 
by comparing said detection amplitude with a threshold value predetermined for decoding. 

[0080] The decoding method thus makes it possible to find the adaptive partitioning effected when the watermarking 
signal is inserted. 

[0081 ] According to advantageous characteristic of this embodiment in the invention, the threshold value predetermined 
for decoding is less than said predetermined threshold value used during the method of determining a partitioning. 
[0082] This is because, at decoding, the detection test is performed on coefficients which have a priori been modulated 
and possibly made noisy. These coefficients are different from the original coefficients, and the result of the detectability 
test at decoding will be different from that performed in order to determine the partitioning. By choosing different thresholds 
on decoding and coding, it is possible to take into account the tendency for variation in detection amplitude due to an 
increased variance in the coefficients which have been modulated and possibly made noisy. 

[0083] Correlatively, the present invention also concerns a device for inserting a watermarking signal in a digital image 
by modulating coefficients of distinct regions, comprising means for partitioning into distinct regions the digital image 
wherein partitioning means are adapted to effect an adaptative partitioning as a function of a criterion of statistical 
detectability of an information bit inserted on each region. 
[0084] According to a preferred embodiment, this device includes: 

means for generating a centred pseudo-random sequence equal in size to the cardinal number of the said set of 
coefficients, formed by centred pseudo-random sub-sequences; 

means for modulation of a set of coefficients by the said centred pseudo-random sequence in order to insert the 
same information bit on the said set of coefficients; and 

means for checking, over each region, a criterion of detectability of the information bit inserted by modulation. 

[0085] Correlatively, the present invention also concerns a decoding device adapted to implement the decoding method 
according to the invention. 

[0086] The present invention also concerns a computer adapted to implement the method of inserting a watermark 
and the decoding method in accordance with the invention. 

[0087] This computer has characteristics and advantages similar to those described with reference to the methods 
which it implements. 

[0088] The present invention also relates to an information storage means which can be read by a microprocessor, 
possibly removable, which stores a program implementing the method of inserting a watermark and/or decoding in 
accordance with the invention. 

[0089] Other particularities and advantages of the invention will also emerge from the following description of the 
invention. 

[0090] In the accompanying drawings, given by way of non-limitative examples: 

Figure 1 is an overall diagram illustrating a device for inserting a watermark into a digital image; 
Figure 2 schematically illustrates a step of spectral decomposition employed when inserting a watermark; 
Figure 3 is an algorithm of the step of determining a partitioning in accordance with afirstembodimentof the invention; 
Figure 4 is an example of partitioning of a region by division; 

Figure 5 is an algorithm of the step of determining a partitioning in accordance with a second embodiment of the 
invention; 

Figure 6 is an algorithm of the decoding method according to a first embodiment of the invention; 

Figure 7 is an algorithm illustrating the application of the method of inserting a watermark of Figure 1 toadigital image; 

Figure 8 is an algorithm illustrating the application of the decoding method of Figure 4 to a digital image; 

Figure 9 is a block diagram illustrating a variant embodiment of the step of determining a partitioning in accordance 

with a third embodiment of the invention; 

Figure 1 0 illustrates the application of the step of determining a partition to a digital image according to a preferred 
embodiment of the invention; 

Figure 11 is an algorithm of the step of determining a partition in accordance with a first way of the preferred 
embodiment of the invention; 

Figure 12 is an algorithm of the step of determining a partition in accordance with a second way of the preferred 
embodiment of the invention ; 

Figure 13 is an algorithm of the insertion method according to one embodiment of the invention; 
Figure 14 illustrates, by means of a graph, the insertion method of the figure 13; and 

Figure 15 is a block diagram illustrating a device adapted to implement the methods in accordance with the invention. 
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[0091] A description will first of all be given, by reference to Figure 1, of the insertion of a watermark into a set of 
coefficients representative of a digital image I. 

[0092] This digital image I can be represented by a set of coefficients, either in the spatial domain or in a transformed 
domain, the coefficients in this latter case being hybrid coefficients, that is to say coefficients located both in the spatial 
5 domain and in the frequency domain. Such a representation of the image is obtained, for example, by using a sub-band 
originating from a space-frequency decomposition of the image, for example a discrete wavelet decomposition. 
[0093] Here, the insertion of the watermark is achieved via a spread-spectrum insertion technique, by modulation of 
coefficients of a space-frequency representation of the image, obtained by the space-frequency transformation of the 
digital image I. 

10 [0094] By way of example, a discrete wavelet decomposition S can be used, represented diagrammatically in Figure 
2. This space-frequency decomposition is well known in the image-processing field, and only the principle thereof is 
reiterated below. It makes it possible to segment the image into sub-bands of frequencies and to obtain hybrid coefficients, 
that is to say spectral coefficients which are also located in space, here in the plane of the image. 
[0095] The image I consists of a series of digital samples. The image I is represented, for example, by a series of 
is bytes, each byte value representing one pixel of the image I, which may be a black-and-white image with 256 grey levels. 
[0096] The multi-resolution spectral decomposition means consist of a sub-band decomposition circuit, or analysis 
circuit, formed by a set of analysis filters, associated respectively with dividers-by-two. This decomposition circuit filters 
the image signal I in two directions, into sub bands of high spatial frequencies and of low frequencies. The circuit includes 
several successive analysis blocks for breaking down the image I into sub-bands according to several levels of resolution. 
20 [0097] By way of example, the image I is broken down here into sub-bands with a level of decomposition d equal to 3. 
[0098] A first analysis block receives the image signal I and filters it through two respectively low-pass and high-pass 
digital filters, in a first direction, for example horizontal. After passing through dividers-by-two, the resultant filtered signals 
are in their turn filtered by two respectively low-pass and high-pass filters, in a second direction, for example vertical. 
Each signal is again passed through a divider-by-two. Hence four sub-bands LL 1? LH 1; HL 1 and HH 1 are obtained at 
25 the output of this first analysis block, with the highest resolution in the decomposition. 

[0099] The sub-band LL 1 includes the low-frequency components in the two directions of the image signal I. The sub- 
band LH 1 includes the low-frequency components along a first direction and high-frequency components along a second 
direction of the image signal I. The sub-band HL 1 includes the high-frequency components along the first direction and 
the low-frequency components along the second direction. Finally, the sub-band HH 1 includes the high-frequency corn- 
so ponents along the two directions. 

[0100] A second analysis block in its turn filters the sub-band of low frequencies LL 1 so as, in the same way, to supply 
four sub-bands LL 2 , LH 2 , HL 2 and HH 2 with intermediate level of resolution in the decomposition. A third analysis block 
then filters the sub-band of low frequencies LL 2 so as to supply four sub-bands LL 3 , LH 3 , HL 3 and HH 3 . Finally in this 
example, the sub-band LL 3 is in its turn analysed by a fourth analysis unit in order to supply four sub-bands LL 4 , LH 4 , 
35 HL 4 and HH 4 with the lowest resolution in this decomposition. 

[0101] Hence thirteen sub-bands and four levels of resolution are obtained. Obviously, the number of levels of reso- 
lution, and consequently of sub-bands, may be chosen differently. 

[0102] Obviously, other types of spectral transformation could be used, such as the discrete Fourier transform, the 
discrete cosine transform or the Fourier-Mellin transform. 
40 [01 03] I n a general way, frequency sub-bands are obtained forming a set of spectral coefficients into wh ich a watermark 
can be inserted. 

[0104] The high-frequency sub-band of highest resolution HH 1 can be considered here in order to insert a watermark. 
This sub-band HH 1 thus supplies a set of coefficients X, equal to N in size, for example. 
[0105] Let X = {X h 1 < i < N} be the set of coefficients representative of the digital image I. 
45 [0106] Obviously, it would be possible to consider a subset of coefficients of this set X, representative only of a part 
of the digital image I. 

[0107] In order to insert a watermark into this set of coefficients X, one technique consists in inserting a pseudo- 
random signal, by spreading its spectrum, so as to render this signal undetectable by spectral or statistical analysis. 
[0108] Let w = {Wj, 1 <j < P} be a watermark of size P < N, a pseudo-random signal with known distribution and average 
50 of zero. The most widespread distributions are the binary distribution {-1, 1}, the uniform distribution over [-1, 1] and a 
centred normalised Gaussian distribution N (0, 1). 
[0109] The linear insertion formula is : 

55 X j = Xj + b otj Wj , with 1 < j < P, 

in which {X;, 1 < j < P} is a subset of the set of coefficients X, b is an information bit, and as is a weighting coefficient, 
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also called modulation amplitude. 

[01 10] The set {j, 1 < j < P} is also called information-bit insertion carrier. 

[0111] In order to insert a signal which may be formed from several information bits, it is worth partitioning the set of 
coefficients so as to determine the number of regions or the number of insertion carriers available in this set. 
[0112] To do that, the insertion device 10 in accordance with the invention comprises a device for determining a 
partition 1 1 and insertion means 12. 

[01 1 3] The device for determining a partition 1 1 is suitable for carrying out adaptive partitioning of the set of coefficients 
as a function of the detectability of an inserted watermark. 

[01 1 4] This method of adaptive partitioning consists in varying the number of coefficients modulated in order to insert 
an information bit of the watermark, so as simultaneously to satisfy the criteria which are set of invisibility and of correct 
detection probability. 

[01 15] A description will first of all be given, with reference to Figure 3, of a method of determining a partitioning into 
distinct regions of a digital image which makes it possible to maximise the number of regions obtained in order to 
maximise the number of information bits inserted in the image. 

[0116] Consideration is given to a representation of the digital image by a set of coefficients, either in the spatial 
domain, or in a transformed domain, the coefficients being in the latter case hybrid, that is to say located both in the 
spatial domain and in the frequency domain. Such a representation of the image is for example obtained by using a sub- 
band issuing from a spatio-frequency decomposition of the image, as a discrete wavelet decomposition previously 
described. 

[01 17] The partitioning of the image is effected in order to insert a watermarking signal by modulation of the coefficients 
of each region. 

[01 18] According to the invention, this partitioning into regions is performed by an adaptive partitioning as a function 
of a criterion of detectability of an information bit inserted on this region. 

[0119] The number of regions obtained under constraint of the detectability criterion makes it possible to define a 
maximum capacity of the digital image. This maximum capacity corresponds to the number of information bits which it 
is possible to insert in the digital image whilst complying with the detectability criterion. 

[0120] Initially, at a step E100, this capacity Q is initialised to 0 and the starting point is an initial partitioning which is 
limited here to a single region R of size P. 

[01 21 ] Naturally, the initial partitioning could be different and notably include already several distinct regions with sizes 
which might be different. 

[01 22] It is assumed that the size P of this initial region R is greater than a minimum size denoted P stat which corresponds 
to the minimum size of a statistically significant sample of coefficients for detecting a watermarking signal inserted in 
these coefficients. Typically, this minimum size can be fixed at 100. 

[0123] This initial region R can consist of the set of coefficients of the representation under consideration or a subset 
of this representation. 

[0124] A test step E101 then tests the detectability of the insertion of a watermarking signal in the initial region R. 
[0125] This detectability test can be effected by actual modulation and detection operations on the coefficients of the 
initial region R. 

[0126] In practice, all the coefficients of the initial region R are copied into a working memory so as not to directly 
watermark the image. 

[0127] Next all the coefficients of the initial region R are modulated by inserting the watermarking signal. 
[0128] Conventionally, the following procedure is followed in order to effect this modulation: 

The initial region R consists of a set of coefficients: 

X:{X|, « = 1,...P} 

Normally a linear modulation is used of the type: 



X i * = X i ±a i W ij with i = 1, ...,P 

where W : {Wj, i = 1, ...,P} is the watermarking signal generally chosen as a pseudo-random sequence with a known 
distribution and a nil mean. 

[0129] The value of the sign ± of the modulation depends on the binary value to be inserted: for example the sign - 
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corresponds to the value 0 and the sign + corresponds to the value 1 . 

[01 30] It is possible to choose, by way of example, a pseudo-random sequence W which follows a uniform law on the 
interval [-1, 1]. 

[0131] Naturally, any pseudo-random watermarking with known distribution and nil mean can be suitable. The most 
usual distributions for the watermarking signal W are, apart from the aforementioned uniform distribution [-1, 1], the 
binary distribution {-1,1} and the centred standardised Gaussian distribution N(0, 1 ). 

[0132] The modulation can optionally be protected by means of a secret key characterising for example the kernel for 
reproducing the pseudo-random sequence W. 

[0133] The term oq is the modulation amplitude, or weighting coefficient, applied to the coefficient Xj. 
[0134] A constant weighting coefficient can be used for all the coefficients of the region R, such that <Xj = a v for any i, 
the value of the weighting coefficients oc v being equal to a maximum value ensuring imperceptibility of the watermarking 
signal for the modulation of the set of coefficients representing said digital image. 

[0135] This maximum value Oy (referred to as Just Noticeable Difference or JND amplitude) decreases in fact with 
the length of the sequence of the watermarking signal W, and consequently with the number of modulatable coefficients. 
[0136] Given that the capacity Q will be maximised when all the coefficients of the representation of the digital image 
are used for the modulation, it is judicious to determine a modulation JND amplitude for a number of modulatable 
coefficients equal to the cardinal number of the set of coefficients representing the image. 

[0137] This maximum value oc v or JND amplitude corresponds to the maximum modulation amplitude which can be 
used in the aforementioned linear insertion mode! beyond which an observer is capable of visually detecting a change 
in the reconstructed watermarked image. 

[0138] In order to determine this JND amplitude, it is possible to use a visibility model which makes it possible to 
predict the visibility of a watermarking operation according to different parameters, namely the representation of the 
signal through a spatio-frequency transformation S used, the sub-band considered for the insertion, the type of distribution 
of the sequence of the watermarking signal W and the length P of the sequence W. 

[01 39] A simple model, developed by WATSON and described in the article entitled "Visibility of wavelet quantization 
noise", A.B. WATSON et al, IEEE Trans, on Image Process, 6(8), 1164-1175, 1997, makes it possible to predict the 
visibility of a set of modulated coefficients from the measurement of visibility of a single modulated coefficient. Reference 
can advantageously be made to this document for the detailed description of this model. 

[01 40] Thus use is made of afunction dependenton the length P of the watermarking signal W, the type of transformation 
S used and the (base) sub-band considered for the insertion, but independent of all the coefficients to be watermarked X. 
[0141] This calculation function used can be written: 



where oc base (S) is a base value, dependent on the transformation used S and on the (base) sub-band considered for the 
insertion, of the maximum weighting coefficient ensuring imperceptibility during modulation of a single coefficient of this 
sub-band, 

P is strictly greater than 2, and 

is the mathematical expectation of the function |tv|P. 
[0142] The base values oc base (S) can be measured once and for all, for each sub-band of coefficients in the wavelet 
decomposition, from a single psychovisual measurement, and be stored in a table of visual amplitudes. 
[0143] p is the exponent of a Minkowsky sum and can be chosen equal to 5, for example. 
[0144] The mathematical expectation £[|w|P ] corresponds to an estimation of the mean of the function |w|P. 
[01 45] This visibility model does not take into account the image I to be watermarked itself and is independent thereof. 
This is equivalent to considering that the image I is uniform. It is a question of a "worst case" model since the presence 
of the image signal makes it possible to visually mask the modulation itself. 

[0146] The JND amplitude can thus be calculated for a length P of the watermarking signal W equal to the cardinal 
number of all the modulatable coefficients of the sub-band under consideration. 

[01 47] However, it is judicious to exploit the fact that the image signal itself makes it possible to mask the modulation. 
[0148] Thus, in this example, use is made, for each coefficient Xj to be modulated, of a weighting coefficient in ac- 
cordance with a law of the form a, = kj.oCy, where k, is a modulation factor dependent on the coefficients situated close 
to the coefficient considered on the region R and oc v is equal to the maximum JND value of the weighting coefficient. 
[0149] Thus each coefficient will be modulated according to the local content, which makes it possible to locally 
increase the amplitude of the modulation to the benefit of detectability. 
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[01 50] In order to test detectability in the test step E1 01 , adetectability amplitude is next calculated from the detectability 
criterion and this detectability criterion is validated by comparing the detectability amplitude with a predetermined thresh- 
old value denoted Tc. 

[0151] A calculation of correlation C(X*, W) between all the modulated coefficients X* and the marking signal W is 
generally made and it is decided whether there has actually been watermarking of the image when the result of the 
correlation calculation C(X*, W) is greater than the predetermined threshold value Tc. 

[01 52] A standardised statistical test can also be used for the detection, such as the one described in the article entitled 
"A method for signature casting on digital images" by I. PITAS, in Proc. ICIP, pages 215-218, September 1996. The 
detection is then characterised in terms of probability. It is thus possible to choose a threshold value Tc corresponding 
to a fixed detection probability level, for example 99.95%. 

[0153] Thus, for a given region R defining a set of coefficients X of size P and with a fixed weighting coefficient a, it 
is possible to calculate a detectability amplitude and compare it with a threshold value Tc so as to validate or not the 
detectability criterion on this region R. 

[0154] The detectability criterion is in fact in theory a function depending on the size P of the region R, the detection 
threshold Tc related to the required detection probability, the watermarking signal W, the variance a 2 x of the coefficient 
X to be modulated on the region R and the weighting coefficient used a for modulating the coefficients X according to 
the linear model described above. 

[0155] Thus, for example, for a detection by standardised hypothesis test and a value of the weighting coefficient a 
which is constant and equal to oc v for all the coefficients Xj of the region R, a minimum length P min (X) is defined by the 
function: 



where a, b, c and d can be assimilated to constants. The signal being non-stationary, that is to say the variance a 2 x 
being a function of space, it is indeed necessary to adapt the modulation length P generated by the watermarking signal 
to the image signal itself. 

[0156] On decoding, the set of coefficients X being modulated, this function also depends on the variance o 2 w of the 
inserted watermarking signal W. 

[0157] In addition, it is necessary to provide for any distortions which the image may undergo after modulation, due 
for example to a compression of the image. 

[0158] When this noise can be considered to be additive and decorrelated from the coefficients, it is possible to use 
a model dependent on the variance o 2 q of this additive noise modelling any distortions applied to said coefficients after 



[0159] Thus, to a given weighting coefficient a, making it possible to ensure invisibility of the inserted watermarking 
signal, a detection threshold Td fixed as a function of the detection probability required at decoding and a set of selected 
coefficients X, there also corresponds a minimum modulation length P min (X) which makes it possible to ensure detecta- 
bility of the inserted watermarking signal. 

[0160] On decoding, for a detection by standardised hypothesis test, this minimum length P m j n (X) is then defined by 
the function: 



where a, b, c and d can be assimilated to constants. 

[0161 ] It is therefore necessary, if the same detectability test is used during coding and decoding, to choose athreshold 
value Tc on coding greater than that Td used at the time of decoding if it is wished to be able to find the same partitioning 
on decoding of the image. 

[0162] In practice, the threshold value Td on decoding which it is wished to use is fixed, that is to say for example the 
required detection probability level when the detection is made by a standardised hypothesis test, and a higher threshold 
Tc on coding is chosen, in order to take account of the variance in the watermarking signal W and the variance of an 
additive noise fixed a priori so as to anticipate the influence of the subsequent processing undergone by the image. 
[0163] This condition thus makes it possible to find the same detectability criterion which constrains the partitioning 
of the digital image without having to transmit to the decoder the partitioning effected at the time of insertion of the 




modulation. 



Pmin(X) = 




\b(<r 2 x +or 2 ov +<y*) + ca 2 +d\ 
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watermarking signal W. 

[0164] In this example embodiment, a partitioning by iterative division of the digital image is effected. It is consequently 
considered, at step E101, that the detectability criterion is satisfied when the calculated detection amplitude is greater 
than or equal to the predetermined threshold Tc. 
5 [0165] If the detectability criterion is not satisfied, then it is probable that the detectability criterion is also not satisfied 
in subregions of the initial region R of smaller size because the detection probability generally decreases with the size 
of the regions. A step E102 of stopping the method of determining a partitioning makes it possible to stop the process 
of iterative division of the image. 

[0166] Insertion not being judged possible in the initial region R, the capacity Q of the image is therefore nil. 
to [0167] If on the other hand the detectability criterion is satisfied, that is to say here that the detection amplitude is 
greater than or equal to the threshold value Tc, it is possible that it is also satisfied for subregions of the initial region R. 
[0168] It is therefore judicious to attempt to divide the region R as long as the detectability criterion is satisfied in order 
to increase the capacity Q. 

[0169] There are of course several possible means for effecting the spatial division of the region R into subregions 
15 Rj, of sizes Pj, with j = 1, K. 

[0170] In this example a division into a quadtree is chosen, as illustrated in Figure 4. 
[0171] This division consists of dividing each region into four subregions of identical size. 

[0172] Thus, by way of example, in Figure 2, a digital image of initial size N is being partitioned into two sub-regions 
R5 and R6 of sizes N/4, seven subregions R1, R2, R3, R4, R7, R8, R9 of sizes N/16 and four subregions R10, 
20 R 11 , R 1 2, R 1 3ofsizesN/64. 

[0173] Thus, at the division step E103, the region R is divided into subregions Rj, j = 1,..., K of sizes Pj. In a case of 
a division into a quadtree, K = 4 and Pj = P/4. 

[01 74] At step E1 04 a first subregion of index j (for example j = 1 ) is considered and a value of potential gain in capacity 
denoted G Q is initialised to -1 . 

25 [0175] The partitioning determination method then includes a prior step E105 of comparing the size Pj of the region 
Rj with the minimum size P stat corresponding to the minimum size of a statistically significant sample for the detection 
of a watermarking signal inserted in said region as described above. 

[01 76] If the size Pj of the region Rj is less than this minimum size P stat , at step E1 08 the index J = J+1 is incremented 
in order to process the following sub-region R j+1 . 
30 [0177] Otherwise, the detectability criterion is checked on the subregion Rj in the same way as previously described 
with reference to step E1 01 . 

[0178] In practice, the watermarking signal is inserted on the coefficients of the said region Rj, a detection amplitude 
is calculated and this detectability amplitude is compared with the threshold value Tc. 

[0179] The detection amplitude can be a correlation value between the modulation signal X* and the carrier of the 
35 watermarking signal W or the result of a standardised hypothesis test. 

[01 80] If the detectability criterion is verified, the value of potential gain in capacity G Q is increased by one unit at step 
E107. 

[0181] Then at step E108 the index j of the subregions is incremented and it is checked at the test step E1 09 that this 
index j is less than or equal to the number of K. 
40 [0182] In the affirmative, steps E105 to E109 are reiterated on the following subregion. 

[0183] Otherwise, when the index j = K+1, it is checked in a test step E1 10, whether the potential gain in capacity 
value G Q is strictly positive, that is to say whether there exist at least two subregions Rj of the initial region R for which 
the detectability criterion is validated. 

[0184] In the negative, the partitioning of higher rank carried out on the initial region R is not validated and this region 
45 r is preserved. Such is the case when the detectability criterion has not been validated for any of the subregions Rj or 
only for a single subregion Rj. 

[0185] In the latter case, it is preferable to keep the region R for inserting the watermarking signal since it is of size P 
greater than the size of the single validated subregion Rj. 

[0186] At step E112, the value of the capacity Q is then increased by one unit and the region R of size P in the 
50 partitioning is kept. 

[0187] An actual insertion step E1 13 by modulation of the coefficients of the region R can then be implemented in a 
normal fashion using the linear insertion model described above. 

[01 88] An elimination step E1 1 4 then makes it possible to eliminate the region R from astack of regions to be processed. 
[0189] On the other hand, if at step E110 the value of the potential gain capacity Gc is strictly positive, in an addition 
55 step E1 1 1 all the subregions Rj for which the detectability criterion is validated in the stack of regions R to be processed 
are added. 

[0190] A return step E 1 15 makes it possible to reiterate all the steps E103 to E1 13 in order to process the regions of 
the stack in a recursive fashion. 
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[0191] At the end of the recursive partitioning method just described, a capacity value Qm corresponding to the 
maximum capacity of the image for insertion of a watermarking signal is obtained, guaranteeing a correct detection of 
this signal. This maximum capacity Qm is equal to the cardinal number of the distinct regions obtained by the said 
partitioning on which the detectability criterion is validated. 
5 [0192] In addition, the insertion step proper E1 13 can be omitted during a first "pass" which then makes it possible 
only to determine the maximum capacity Qm of the image and the validated regions. 

[0193] Once the maximum capacity Qm has been calculated and the corresponding partitioning is stored in memory, 
it is then possible for a user to determine the message which he wishes to insert, corresponding to a number of information 
bits equal to the maximum capacity Qm, and then to insert this message in a normal manner by modulating the coefficients 
10 of the different regions obtained by the previously determined partition. 

[0194] A description will now be given, with reference to Figure 5, of a second embodiment of the invention in which 
it is wished to insert a message of fixed capacity Qs in an image I. Naturally, this fixed capacity Qs must have a value 
less than the maximum capacity Qm of the digital image to be watermarked. 

[0195] As illustrated in Figure 5, the method of determining a partitioning in this embodiment includes steps E100 to 
15 E113 identical to those of the method of determining a partitioning of the embodiment of Figure 3 and which it is not 
necessary to describe again here. 

[0196] As before, at step E1 12, the value of the capacity Q is incremented when the partitioning of rank directly higher 
is not validated. 

[01 97] Then the method includes a comparison step E1 1 6 in which this value of the capacity Q is compared with the 
20 fixed capacity value Qs. 

[0198] In practice, it is tested at step E1 1 6 whether the capacity value Q is equal to the fixed capacity value Qs. 
[0199] In the affirmative, a stop step E1 17 makes it possible to stop the iterative partitioning of the image since there 
are a number of regions equal to the number of information bits Qs which it is wished to insert. 

[0200] Otherwise, in an elimination step E1 14, the processed region is eliminated from the stack of regions R to be 
25 processed formed at step E1 1 1 . 

[0201] The method of determining a partitioning also includes, in this embodiment, a sorting step E1 18 which sorts 

the regions R to be processed in the stack by detectability amplitude value calculated on these regions R. 

[0202] Thus a return step E1 1 5 reiterates the partitioning from the division step E1 03 for a region R in the stack having 

the highest detectability amplitude amongst all the other regions to be processed. 
30 [0203] A sorting variant can consist of sorting the regions to be processed in decreasing order of size so as to process, 

as a priority, the largest regions which a priori provide better detectability. 

[0204] Naturally, a partitioning by iterative merging of the digital image could be carried out in the same way, for 
example also using a quadtree partitioning structure. 

[0205] The initial partitioning then consists of regions with smaller acceptable sizes, greater than or equal to the 
35 previously defined minimum size P stat . 

[0206] For each region of the partitioning for which the delegability criterion is not validated, this region is merged 
with another region of the partitioning. 

[0207] This other region of the partitioning is if possible a region for which the detectability criterion is not validated 
so as to increase the capacity. 

40 [0208] Alternatively, another region can be a region satisfying the detectability criterion so as to make detection reliable. 
[0209] There are many possibilities for the choice of the region to be merged. It is possible to take as the merger 
criterion the detection amplitude calculated on each of the regions and merge for example the regions associated with 
the lowest detection amplitudes. 

[0210] In this particular case of a merging in a quadtree partitioning structure, the merging is structurally constrained 
^5 and consists of deciding whether or not four adjacent blocks, for example R1 , R2, R3 and R4 in Figure 4, must be merged. 
[0211] A description will now be given, with reference to Figure 6, of the associated decoding method which makes 
it possible to find the partitioning when the watermarking signal inserted in a digital image is extracted. 
[0212] Indeed, because of the adaptation of the partitioning to the digital image, this partitioning will not be known a 
priori at the decoder. 

50 [0213] The decoding method then consists, in general terms, of testing all the acceptable partitionings. The same 
iterative implementation technique is used as the one used when the partitioning is determined for inserting the water- 
marking signal. 

[021 4] Thus a description will be given of an example embodiment in the case of a segmentation by dividing by blocks 
in aquadtree, symmetrical to that using the partitioning determination method of the embodimentdescribed with reference 
55 to Figure 3. 

[0215] At an initialisation step E200, there is an initial partitioning limited to a single region R of size P. This region 
consists of all the coefficients of the same representation of the digital image as the one used at the time of insertion of 
the watermarking signal. 



12 



EP 1 043 687 B1 



[0216] A test step E201 makes it possible to verify that the size P of the initial region R is indeed greater than the 
minimum size P stat as described previously. 

[0217] In the affirmative, at step E203 a criterion of detectability of the information bit inserted on the region R, obtained 
by demodulation of the coefficients on this region R, is verified. 
5 [0218] In practice, a detection amplitude is calculated for all the coefficients of the region R and the detectability 
criterion is validated by comparing this detection amplitude with a predetermined threshold value Td for decoding. 
[0219] It is possible to proceed by means of a standardised statistical test, and the result compared with the prede- 
termined threshold value Td. 

[0220] As explained before, it is necessary, if the same detectability test is used during coding and decoding, to choose 
to a threshold value Tc on coding greater than that Td used at decoding if it is wished to be able to find the same partitioning 
on decoding of the image even when the latter has been made noisy by post-processing. 

[0221] This condition thus makes it possible to find the same detectability criterion which constrains the partitioning 
of the digital image without having to transmit, to the decoder, the partitioning effected when the watermarking signal W 
was inserted. 

15 [0222] If the detectability criterion is not verified at step E203, a division is carried out, at step E204, of the region R 
into subregions Rj and sizes Pj with j = 1 , K. Then the same partitioning structure is used as the one used on coding, 
that is to say, in this example, a division of the region R into four subregions Rj of identical sizes. 
[0223] Then, in an addition step E205, the subregions Rj are added in a stack of regions R to be processed and, at a 
return step E202, the decoding method is reiterated recursively as from the test step E201 on each region R of the stack. 

20 [0224] When at the step of calculating the detectability criterion E203, the latter is verified for a region R, there is 
extracted at the extraction step E206 the information bit inserted in this region R as an element of the method inserted 
in the digital image I. This extraction is effected in a conventional fashion by calculating a detection amplitude. The result 
of the test gives, in absolute value, a detection amplitude to be compared with the threshold value Td of the decoding, 
the sign giving the value 0 or 1 of the information bit inserted. 

25 [0225] It should be noted that this extraction is already partly effected during the detection step E203 so that the 
detection E203 and extraction E206 steps are structurally coupled. 

[0226] An elimination step E207 is implemented in order to eliminate the region R from the stack of regions to be 
processed, which was formed at the addition step E205, and the return step E202 then makes it possible to implement 
the decoding method recursively on each region R of the stack to be processed. 
30 [0227] A description will now be given, with reference to Figures 7 and 8, of a practical application of the method of 
determining a partitioning of a digital image in a particular embodiment of the invention. 

[0228] In this application, the watermarking signal is inserted by means of a technique of insertion by spectrum spread- 
ing, by modulation of coefficients of a space-frequency representation of the image, obtained by a spatio-frequency 
transformation S of the digital image I as previously described with reference to figure 2. 
35 [0229] In general terms Sj sub-bands are obtained, with J = 1 , M and M equal hereto 13 at a transformation step E300. 
[0230] Next, for each sub-band Sj of index j considered at step E301, a partitioning is effected at step E302 by 
implementing the method of determining a partitioning described previously to with reference to Figure 3. 
[0231] In this way a capacity Qj is derived for each sub-band Sj. 

[0232] A step E303 of insertion by modulation of the coefficients next makes it possible to insert Qj information bits in 
^o the Qj regions validated. 

[0233] In steps E304 and 305, the following sub-bands are next considered as long as all the sub-bands Sj have not 
been processed, and the steps E302 and E303 of partitioning and inserting a watermarking signal of capacity Qj are 
reiterated on each sub-band Sj. 

[0234] Next, an inverse transformation S _1 is applied to the image in order to obtain the watermarked image, of total 
45 capacity Q equal to the sum of the capacities Qj calculated for each of the sub-bands Sj. 

[0235] Thus, each sub-band of the representation is considered independently for the adaptive partitioning which is 
carried out both at coding and at decoding. 

[0236] As illustrated in Figure 8, the decoding method in this embodiment is implemented in a similar fashion using 
the decoding method described with reference to Figure 6, the initial region R here corresponding each time to a sub- 
50 band of the decomposition of the image. 

[0237] A step E400, identical to step E300, divides the image I into Sj sub-bands with j = 1 , ... M. 

[0238] A first sub-band Sj is considered at step E401 . The order of processing of the sub-bands is identical to that 

used during coding in order to extract the inserted message in order. 

[0239] The adaptive partitioning is implemented at the partitioning step E402 in order to calculate the capacity Qj of 
55 the sub-band Sj. Next, at the extraction step E403, the Qj information bits are extracted by demodulating the coefficients 
of the Qj regions of the partitioning. 

[0240] At the test steps E404 and E405, it is checked in the usual fashion that all the sub-bands have been processed 
and, in the negative, the decoding steps E402 and E403 (partitioning and extraction) are reiterated on the remaining 
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sub-bands. 

[0241] In this way, at the output of the decoder, a message of Q bits is obtained, Q being equal to the sum of the 
capacities Qj of each sub-band Sj. 

[0242] In a variant of this method illustrated in Figure 9, it is possible to effect a partitioning on a post-processed image 
signal in order to obtain an a priori robustness of the watermarking signal to any post-processing undergone by the image. 
[0243] The principle of this variant consists of performing a detectability test in a test step E309 which is no longer 
coupled to the steps of partitioning E302 and insertion E303 described previously. 

[0244] After the insertion of the watermarking signal in a region, a reverse transformation S _1 is first of all performed 
at step E306 in order to find the watermarked image. Then, in a step E307, post-processing is applied, such as a 
compression of the digital image. 

[0245] Next, in a step E308, the transformation S is reiterated in order to find the watermarked coefficients, possibly 
made noisy by the post-processing carried out on the image at step E307, and the detectability criterion is calculated 
as previously described. 

[0246] Steps E300 to E309 are reiterated on another acceptable partitioning of the image, here for example by dividing 
the relevant region R into four, and the detectability test is performed at step E309 on the four sub-regions. 
[0247] This partitioning is adopted only if it represents a gain in capacity, that is to say, in this example of quadtree 
division, if the detectability criterion is validated on at least two sub-regions. 

[0248] This variant embodiment makes it possible to directly take into account the post-processings which will be 
applied to the digital image. 

[0249] A method of determining a partition according to a preferred embodiment of the invention will now be described, 
by reference to Figure 10. 

[0250] The method of inserting a watermark is employed on a digital image I. 

[0251] This insertion method first of all includes a method of determining a partition of a set of coefficients representative 
of the digital image I. 

[0252] This partition is applied as a function of a criterion of detectability of an information bit inserted on each region 
of the partition. 

[0253] It should be noted that, in practice, a partition of a digital image into a maximum number of regions satisfying 
a detectability criterion can be sought so as to ascertain the practical capacity of this image with regard to the insertion 
of a watermark. This maximum capacity of the digital image corresponds to the 

[0254] number of information bits which it is possible to insert into the digital image while complying with the detectability 
criterion. 

[0255] It is also possible, for a given watermark, comprising a given number of information bits of less than the maximum 
capacity of the image I, to seek the best possible partition of the image I which makes it possible to insert this watermark 
of given length, for a given detection probability, and in regions of the largest possible sizes so as best to ensure invisibility 
of the inserted watermark. 

[0256] The method of determining a partition first of all includes a step of spectral decomposition E500 of the wavelet- 
decomposition type described previously (Discrete Wavelet Transform or DWT). 

[0257] Thus a set of spectral coefficients X is obtained, representative of the digital image I . These spectral coefficients 
are distributed into frequency sub-bands as illustrated diagrammatically in Figure 2, on completion of the decomposition 
step E500. The size of the set of coefficients X is equal to N = M x M for a square image. Obviously, the same method 
could be applied to a rectangular image. 

[0258] In this embodiment of the invention, the method includes a step E501 of generating a centred pseudo-random 
sequence of size N equal to the cardinal number N of the set of coefficients X representative of the image I. 
[0259] Let w = {w mn , 0 < m, n < M} be this pseudo-random sequence. 

[0260] This pseudo-random sequence w is formed from centred pseudo-random sub-sequences w k . 

[0261] The pseudo-random sequence can be seen as the merging of the pseudo-random sub-sequences w k . 

[0262] Here, the centred pseudo-random sequence w is formed from pseudo-random sub-sequences w k of identical 

size which correspond to a block, called unit block, of the set of coefficients X. 

[0263] Obviously, the pseudo-random sub-sequences w k could differ in size from one another. 

[0264] Here, the size of these unit blocks corresponds to the minimum number of coefficients suitable for forming a 

region of a size which is statistically significant for satisfying the detectability criterion over this region. 

[0265] As described later, when a standardised statistical detection test is applied to a set of coefficients, a minimum 

number of coefficients is necessary for this test to be significant. 

[0266] The pseudo-random sub-sequences w k are of a size A = L x L which is greater than or equal to this minimum 
number of coefficients. 

[0267] For preference, pseudo-random sub-sequences w k are chosen to be equal in size to this minimum size. This 
characteristic is particularly advantageous when partition is achieved by iterative combination by combining blocks from 
an initial segmentation of the image into unit blocks of size L x L. The initial segmentation then corresponds to the 
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maximum partition of the image which can be applied, beyond which the blocks of coefficients generated are no longer 
statistically significant for the insertion and detection of an information bit. 
[0268] For example, L = 8 is chosen. 

[0269] The pseudo-random sequence w is then the merging of N/64 pseudo-random sub-sequences w k forming unit 
blocks. 

[0270] On each unit block, in accordance with the invention, a centred pseudo-random signal of predefined distribution 
(Gaussian, uniform, etc.) has to be generated on the basis of a function depending on a secret key K and on the index 
k of each unit block. This key can be denoted f(K, k). It is then necessary to ensure that, in practice, the average of the 
signal w k is exactly equal to zero over each unit block. 

[0271] A first method for generating such a centred pseudo-random sequence consists in centring the sequence in a 
deterministic way, by symmetrisation. This method is suitable for generating a pseudo-random sub-sequence of even 
size A. 

[0272] In this case, each centred pseudo-random sub-sequence w k is created by generating one half, A/2, of the 
pseudo-random samples of the sub-sequence via a pseudo-random number generator of known centred distribution, 
and by symmetrising the pseudo-random samples thus generated in order to obtain the other half A/2 of the samples 
of the sub-sequence w k . 

[0273] For a single-dimensional sequence, a set of samples of size A/2 {Wj, 1 < j < A/2} is therefore generated, and 
this set is made symmetric in order to obtain the other half of the samples {Wj = -w A _j, A/2 < j < A}. 
[0274] In order to reinforce the random character of the sequence thus obtained, permutation can be carried out on 
the set of samples obtained {Wj, 1 < j < A} by using a secret key. 

[0275] A second method of generating a centred pseudo-random sequence consists in creating each centred pseudo- 
random sub-sequence w k by generating pseudo-random samples via a pseudo-random number generator and by re- 
distributing the sum of these pseudo-random samples generated over each sample. 

[0276] Thus, first of all, by using the secret key f(K, k), a sequence of pseudo-random samples of the same size A as 

the unit block, i.e. {w i5 1 < i < A}, is generated. 

[0277] Next the exact sum of these samples is calculated: 

A 

S k = ]Tw, 

[0278] Then this sum S k is redistributed over the samples so as to obtain a centred pseudo-random sequence, that 
is to say one with a sum exactly equal to zero. 
[0279] Thus w k = {w r S k /A, 1 < i < A} is obtained. 

[0280] The pseudo-random sequence w is thus composed of the merging of the centred sub-sequences w k in such 

a way that the partitioning of a set of coefficients modulated by the centred pseudo-random sequence w gives rise to 

subsets of coefficients also modulated by centred pseudo-random sub-sequences. 

[0281] Obviously, other techniques could be used to generate centred pseudo-random sequences. 

[0282] As illustrated in Figure 1 0, a step E502 of modulation of the set of coefficients by the centred pseudo-random 

sequence w is then implemented so as to insert an information bit onto the set of coefficients originating from the spectral 

decomposition. 

[0283] In practice, the set of coefficients is copied into a working memory, for example a random-access memory in 
a computer, so as not to mark the image I directly. 

[0284] The same information bit is modulated onto the set of coefficients, for example b = 1 here, according to the 
modulation formula: 



X'i = Xj + oti Wj, with 0 < i < N. 

[0285] The weighting amplitude a, is chosen in the usual way so as to guarantee invisibility of the inserted information bit. 
[0286] A constant weighting coefficient can be used for all the coefficients, such that a, = a v for all i, the value of the 
weighting coefficient being equal to a maximum value ensuring imperceptibility of the watermark for modulation of 
the set of coefficients representative of the digital image. 

[0287] It may be worthwhile exploiting the fact that the image signal itself makes it possible to mask the modulation. 
[0288] To this end, it is possible, for each coefficient Xj to be modulated, to use a weighting coefficient according to 
a rule of the form a, = k,.a v , where k, is a modulation factor depending on the coefficients situated in the vicinity of the 
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coefficient in question and oc v is equal to the maximum value of the weighting coefficient ensuring imperceptibility of the 
watermark for the modulation of the set of coefficients. 

[0289] Thus, each coefficient will be modulated as a function of the local content of the image, which makes it possible 
locally to augment the amplitude of the modulation to the benefit of the detectability. 

[0290] For preference, in this embodiment of the invention, it is wished to apply a partitioning to a post-processed 
image signal so as to obtain robustness, a priori of the inserted watermark, against any post-processing undergone by 
the image I. 

[0291] To do that, a distortion is applied, in a distortion step E503, to the set of modulated coefficients. 
[0292] By virtue of the modulation of the coefficients in a single operation, on the basis of a centred pseudo-random 
sequence, the distortion step E503 can be implemented once only upon determination of a partition, in contrast to the 
prior techniques in which each sub-region of coefficients was modulated independently of the other and required distortion 
to be applied to each level of the partition. 

[0293] This is all the more advantageous since, here, the step E503 of applying a distortion includes the following 
sub-steps as illustrated in Figure 10: 

inverse spectral transformation E503a of the set of coefficients X'j modulated in order to obtain a watermarked image; 
application proper, E503b, of distortion to the watermarked image; and 

spectral transformation E503c of the watermarked image in order to obtain a set of modulated coefficients after 
distortion. 

[0294] The distortions applied are distortions which alter the value of the pixels without altering the geometry of the 
image, such as the addition of noise, alteration of the contrast, compression, low-pass or high-pass filtering, for example. 
[0295] Compression of the image is considered here, according to the JPEG standard (for Joint Photographic Expert 
Group) with specified quality factor, with q = 75 typically by default for JPEG compression. 

[0296] In practice, an inverse wavelet transformation DWT~ 1 is applied here, then compression followed by decom- 
pression of the image, and finally forward wavelet transformation again in order to be able to estimate the detectability 
of the information bit inserted in the spectral domain. 

[0297] Next an iterative partitioning of the set of coefficients is applied, so as to search for a partition of the set of 
coefficients in which each inserted information bit can be detected correctly according to a fixed detection probability. 
The number of regions thus determined during the partition corresponds to the practical insertion capacity of the image I. 
[0298] The means for displaying to the users Q bit determined by calculating the capacity should exist. 
[0299] Various techniques for iterative segmentation by division/combination can be used as previously described. 
[0300] A first way of performing the preferred embodiment of the invention will firstly be described, by reference to 
Figure 1 1 , in which the partitioning is achieved by iterative division on the basis of an initial segmentation of the set of 
coefficients into regions of a statistically significant size for validating a detectability criterion over these regions. 
[0301] Following on from the steps of transformation E500, generation E501, modulation E502 and distortion E503 
as described above, a first frequency sub-band will be considered, in a step E504, as a set of coefficients to be partitioned. 
[0302] The initial segmentation of this set of coefficients corresponds to a spatial segmentation of the image I into 
blocks of unit size as defined above. 

[0303] By its very construction, the pseudo-random modulation sequence is centred on each of these unit blocks. 
[0304] In a quadtree-mode combination, it is always four blocks of equal size which are grouped together, called "child 
blocks" in order to form a "parent block", and this combination is recommenced iteratively, giving the notion of recur- 
siveness. 

[0305] Obviously, partition takes place frequency sub-band by frequency sub-band of the signal arising from the 
spectral decomposition. 

[0306] Each unit block is considered here, in a step E505, in an order of scanning of the sub-band, for example from 
left to right and from bottom to top. 

[0307] A validity criterion is applied over a first unit block which has not yet been processed. 

[0308] In a calculation step E506, a statistical detection test is applied to this unit block of modulated coefficients X'j. 
[0309] An example of a standardised statistical test is given below. 

[0310] Let X'j, 1 < i < P be the set of coefficients to be tested after insertion of the pseudo-random sequence w and 
application of pre-distortion. 

[031 1 ] The test consists in calculating, for this set of size P: 
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in which M represents an average value of correlation between the values of each coefficient X'j and of modulation w h 
and V the variance of these values. 
[0312] Hence in practice: 



[0313] The calculation of this test value T requires knowledge of the pseudo-random modulation sequence w which 
is generated once and for all at the above-described generation step E501 . 

[031 4] The detectability is given by the comparison of this test value T with a certain threshold Pp, which corresponds 
theoretically to a probability of correct detection of p% if the test T follows a Gaussian law with variance equal to 1 for 
the two hypotheses tested (presence or absence of watermark). 

[0315] A comparison step E507 makes it possible to compare this test value T with the threshold value P p or correct 



[0316] On completion of this comparison, if the calculated value T is higher than the correct detection rate, the block 
tested is held in memory, in a memory-storage step E509, as being a region available for inserting an information bit. 
[0317] In practice, a capacity-calculating step E51 0 can be implemented in such a way that, for every block of coef- 
ficients tested which is suitable for having an information bit inserted, the value of the capacity is increased by one unit. 
[0318] In contrast, if the calculated value T is less than the correct detection rate Pp, a block of higher level is selected, 
in a combining step E508, that is to say, here, a parent block corresponding to the combination of four unit blocks, and 
the steps E506 of calculating the test value T and of comparison E507 with the correct detection rate P p are reiterated 
over this new block. 

[0319] During this combining of child blocks in order to form a parent block, those child blocks are preferably chosen 
for which the detectability criterion is not satisfied. 

[0320] For application of this detectability criterion, any standardised statistical test can be used for the detection, 
such as that described in the article entitled "A method for signature casting on digital images" by I. PITAS, in Proc. ICIP, 
pages 215-218, September 1996, which makes it possible to characterise the detection in terms of probability. 
[0321 ] A threshold value P p can be chosen, corresponding to a fixed level of detection probability, 99.95% for example. 
[0322] Hence, for a given region defining a set of coefficients X of size P and with a fixed weighting coefficient a, it is 
possible to calculate a detectability amplitude and to compare the latter with a threshold value P p so as to validate or 
not validate the detectability criterion over this region. 

[0323] In a general way, the detection probability decreases with the size of the regions tested, such that, by combining 
regions in order to increase the size of the block of coefficients used for inserting a watermark, the value of the test T is 
increased over this new block. 

[0324] Hence, on completion of these steps of adaptive partitioning E504 to E508, an image insertion capacity can 
be determined, equal to Q bits, as can an appropriate partition of the image corresponding to this capacity Q. 
[0325] A watermark can then be formed, with length equal to Q bits, and can then be inserted, bit by bit, into the 
various regions of the partition. 

[0326] Obviously, the step E520 of insertion proper of the watermark can be coupled with the determination of the 
partition, in such a way that each region available on conclusion of the detectability test E506, E507, can actually be 
watermarked by the insertion of an information bit by modulation of the coefficients of this region. 
[0327] The insertion step E520 makes it possible actually to insert an information bit of the watermark, by choosing 
a coefficient value b = 1 or = -1 in the above-described modulation formula, the alternating of the +/-1 bits conveying 
the content of the message inserted. 

[0328] A centred pseudo-random sequence of known distribution can be used, different for each region of the partition, 
and initialised by a secret key K', as illustrated in Figure 1 . 

[0329] It should be noted that this secret key K' may, if appropriate, be identical to the secret key K used during the 




and 




detection rate. 
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generation of the pseudo-random sequence w for determining the partition. 

[0330] Moreover, a centred pseudo-random sequence w identical to that generated at the generation step E501 could 
be generated. 

[0331] A partition of an image I is shown diagrammatically by way of example in Figure 10, on conclusion of steps 
5 E504 to E508, in which, for example, 25 separate regions are obtained making it possible to insert a 25-bit watermark. 
[0332] Obviously, any other partitioning technique could be used. 

[0333] A second way of performing the preferred embodiment of the invention will now be described, by reference to 
Figure 12, in which the partitioning is achieved by iterative division of the set of coefficients. 

[0334] On conclusion of the steps of transformation E500, generation E501 , modulation E502 and distortion E503 as 
to described above, a first frequency sub-band will be considered, in a step E51 1 , as set of coefficients to be partitioned. 
[0335] A quadtree-mode division will be considered here, in which one block, called "parent block" is always divided 
into four blocks of equal size, called "child blocks". 

[0336] Obviously, partition takes place frequency sub-band by frequency sub-band of signal originating from the 
spectral decomposition. 

15 [0337] Hence, in a step E512, an initial block corresponding to a first frequency sub-band is considered. 
[0338] This initial block is stored in memory in a list L of blocks to be processed. 

[0339] Next, in a step E513, a first block B to be processed, coming from the list of blocks to be processed L, is 
considered. 

[0340] A detectability criterion is applied on this block B in a checking step E514. 
20 [0341] In practice, as previously described, the value T of a statistical detection test on this block of modulated coef- 
ficients is calculated in a calculation step E514a. This step is analogous to the calculation step E506 described by 
reference to Figure 1 1 . 

[0342] The detectability is then given by the comparison of this test value T with a certain threshold P p , which theo- 
retically corresponds to a correct detection probability of p% if the test T follows a Gaussian law with variance equal to 

25 1 for the two hypotheses tested (presence or absence of watermark). 

[0343] A comparison step E51 4b makes it possible to compare this test value T with the threshold value P p or correct 
detection rate. This comparison step E51 4b is analogous to the comparison step E507 described by reference to Figure 1 1 . 
[0344] On completion of this comparison, if the calculated value T is less than the correct detection rate P p , a check 
is made, in a test step E515, as to whether there are blocks remaining to be processed in the list L. 

30 [0345] If not, the recursive partitioning algorithm is terminated. 

[0346] Otherwise the following block B of the list L is considered at step E51 3, and the checking step E51 4 is reiterated. 
[0347] If, on completion of the comparison step E514b, the calculated value of the test T is higher than the correct 
detection rate Pp, detectability is tested over each child block of this parent block B in a checking step E516. 
[0348] In practice, the parent block B is divided into four child blocks and the detectability criterion is checked on each 

35 of these blocks according to a step analogous to the previously described checking step E514. 

[0349] A test E51 7 makes it possible to verify whether the detectability criterion is valid over at least one child block. 
[0350] if yes, a division step E51 8 is implemented in order to actually divide the block B into four child blocks, and an 
addition step E519 allows these child blocks to be added to the list of blocks to be processed L. 

[0351 ] The test and checking steps E51 3 to E51 9 are then reiterated over each of the blocks of the list L until the latter 
^o is empty. 

[0352] In contrast, if, on conclusion of test E517, the detectability criterion is not satisfied on any of the child blocks 
of the parent block B, the tested block B is held in memory in a memory-storage step E509 as a region available for 
inserting an information bit. 

[0353] In practice, a step of updating the capacity E510 can be implemented as described previously, in such a way 
45 that, for each block of coefficients tested which is suitable for insertion of an information bit, the value of the capacity is 
increased by one unit. 

[0354] Hence, on conclusion of these adaptive partitioning steps E51 1 to E519, an image insertion capacity can be 

determined, equal to Q bits, as can an appropriate partition of the image corresponding to this capacity Q. 

[0355] A watermark can then be formed, with length equal to Q bits, and can then be inserted, bit by bit, into the 

50 various regions of the partition. 

[0356] Obviously, as before, the step E520 of insertion proper of the watermark can be coupled with the determination 
of the partition, in such a way that each block B available on conclusion of the detectability test E507, can actually be 
watermarked by the insertion of an information bit by modulation of the coefficients of this block B. 
[0357] The insertion step E520 is carried out in the same way as when partition is achieved by iterative combination. 

55 [0358] Moreover, in one enhanced embodiment of this technique of partition by iterative division, a division of a parent 
block B into child blocks can be adopted on completion of the detectability test E507 only if it represents a gain in capacity, 
that is to say into a quadtree-mode division only if the detectability criterion is satisfied over at least two child blocks. 
Otherwise, it may be preferable to keep the parent block of larger size so as to guarantee better invisibility of the inserted 
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watermark. 

[0359] Obviously, the partitioning techniques described previously are not in any way limiting. 

[0360] Notably, it is possible to use various iterative segmentation techniques by division/merging, to which there 
belong the quadtree block partitioning described previously, but also graph partitioning. In general terms, there are the 
so-called bottom-up methods and the top-down methods depending on whether the iterative partitioning means is merging 
or division. 

[0361 ] Hence it is possible, if appropriate, to combine two techniques of partitioning by division and combination. This 
is the case when an initial partition is chosen, for example, intermediate between the two extremes which are, on the 
one hand, a single region and, on the other hand, a partition into regions of the smallest admissible size, allowing both 
the operations of division and of combination. 

[0362] It is also possible to implement the method of determining a partition in order to maximise the capacity of the 
digital image by using several different segmentation techniques and then choosing the partition which confers the 
greatest capacity Q for the digital image. 

[0363] It is also possible to determine the best partition possible, that is to say the one yielding a set of blocks of 
coefficients of largest possible size, for inserting a fixed-length message, less than the maximum capacity of the digital 
image to be watermarked. 

[0364] In all cases, the present invention makes it possible to generate and to insert a pseudo-random sequence once 
only over a set of coefficients which are available for carrying out the insertion, and, as a corollary, to apply only one 
single pre-distortion to the set of coefficients. 

[0365] It is thus possible rapidly to obtain a partition of the set of coefficients which makes it possible to insert a 
watermark with a fixed correct detection probability for an entire range of distortion lying between "no distortion" and 
"predetermined distortion" applied during the determination of the partition. Moreover, the simulation of maximum dis- 
tortion before insertion is a quasi-deterministic operation which confers an advantage with respect to a theoretical 
statistical calculation. 

[0366] A description will be given now of the insertion method according to the invention in a preferred embodiment. 
[0367] The watermarking signal is inserted after a transformation step E600, in a set of coefficients of the domain 
transformed into sub-bands. 

[0368] As illustrated in Figure 13, in a choosing step E601 , there is for example chosen the high-frequency sub-band 
HH 1; with the first resolution level, corresponding to a high-pass filtering in the horizontal and vertical directions. 
[0369] There is thus a set of spatio-frequency coefficients, with a cardinal number equal to N, denoted for example X 
= {Xj, 1 < i < N}. 

[0370] For an image I of size 512x512, the sub-band HH 1 is of size N = 256x256. 

[0371] Since the watermarking has to be imperceptible and indelible, and therefore difficult to locate and remove by 
pirating, it is chosen to insert a pseudo-random watermarking signal, spreading its spectrum in order to make this signal 
undetectable by spectral or statistical analysis. 

[0372] For example, a pseudo-random sequence w is considered, which follows a uniform law on the interval [-1,1] 
with w = {w h 1 < i < P}, the length P of the sequence being less than or equal to N. 

[0373] Naturally, any pseudo-random watermarking signal, with a known distribution and a nil mean, may be suitable. 
The most usual distributions for the watermarking signal w are, apart from the uniform distribution on [-1 , 1] mentioned 
above, the binary distribution {-1 ,1} and the centred standardised Gaussian distribution N(0, 1). 

[0374] The modulation which it is wished to apply to the coefficients in order to insert the watermarking signal uses a 
linear model of the type: 



X'j = Xj (i ) + awj, with 1 < j < P 

where 

Xj (i ) is a subset of spectral coefficients chosen from the set of coefficients X, and a is a weighting coefficient. 
[0375] Here, the weighting coefficient a is constant for any j. 

[0376] In this example embodiment, it is sought to determine whether the insertion of the watermarking signal w is 
actually possible in the set of coefficients X whilst meeting predetermined criteria of imperceptibility and detectability. 
[0377] According to the invention, the insertion method first of all includes a step E602 of calculating a maximum value 
oc v of the weighting coefficient a as a function of the length P of the watermarking signal w ensuring imperceptibility of 
the watermarking signal w. 

[0378] This maximum value or visual amplitude oc v corresponds to the maximum modulation amplitude which can be 
used in the aforementioned linear insertion model beyond which an observer is capable of visually detecting a change 
at the level of the reconstructed watermarked image. 
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[0379] In this embodiment, a visibility model is used which makes it possible to predict the visibility of a watermarking 
operation as afunction of the different parameters, which are the representation of the signal through the spatio-frequency 
transformation T used, the sub-band HH 1 considered for the insertion, the type of distribution of the sequence of the 
watermarking signal w and the length P of the sequence w. 

[0380] A simple model, developed by WATSON and described in the article entitled "Visibility of wavelet quantization 
noise", A.B. WATSON et al, IEEE Trans, on Image Process, 6(8), 1164-1175, 1997, makes it possible to predict the 
visibility of a set of coefficients modulated from the measurement of visibility of a single modulated coefficient. Reference 
can advantageously be made to this document for a detailed description of this model. 

[0381] During the step E602 of calculating the maximum value oc v , use is thus made of a function dependent on the 
length P of the watermarking signal w and on the type of spectral transformation T used, but independent of the set of 
coefficients to be watermarked X. 
[0382] This calculation function used can be written: 



where oc base (S) is a base value, dependent on the transformation used S and the (base) sub-band considered for the 
insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of this sub-band is 
modulated, 

p is strictly greater than 2, and 

£[| is the mathematical expectation of the function |w|P. 

[0383] The base values oc base (S) can be measured once and for all, for each sub-band of coefficients in the wavelet 

decomposition, from a single psychovisual measurement, and be stored in a table of visual amplitudes. 

[0384] p is the exponent of a Minkowsky sum and can be chosen equal to 5, for example. 

[0385] The mathematical expectation E[| corresponds to an estimation of the mean of the function \ w\$. 

[0386] This visibility model does not take into account the image I to be watermarked itself and is independent thereof. 

This is equivalent to considering that the image I is uniform. It is a question of a "worst case" model since the presence 

of the image signal makes it possible to visually mask the modulation itself. 

[0387] In this example, this calculation step E602 is implemented in order to calculate the maximum value oc v of the 
weighting coefficient a for a length P of the watermarking signal w equal to the cardinal number N of the set of modulatable 
coefficients X. Thus a maximum value denoted oc v (N) is calculated, which ensures imperceptibility of the watermarking 
signal w at the limit of perceptibility. 

[0388] In the example embodiment described here, for a value of 3 = 5, the visual amplitude for the set of coefficients 
of the sub-band HH 1 is given by: 



[0389] Next, in a calculation step E603, the length P of the watermarking signal w is calculated for a minimum value 
oc D of the weighting coefficient a equal to said maximum value oc v (N) calculated and for a predetermined detection 
probability level. 

[0390] This is because, in order to meet a predetermined detectability criterion, the weighting coefficient a must be 
greater than a minimum value a D referred to as the detection amplitude. 

[0391] At the decoder, the received image I*, which corresponds to the watermarked initial image, possibly made 
noisy because of certain distortions applied to the image, is first of all transformed by a transformation identical to that 
used during the insertion of the watermarking signal in order to find once again the set of coefficients which are modulated. 
[0392] Thus asetX* = {Xj*, 1 < i< N} is obtained, which represents the set of coefficients liable to have been modulated. 
By hypothesis, it is assumed that any distortion suffered by these coefficients can be modelled by an additive noise, 
decorrelated from the image signal itself n = {n j5 1 < i < N}. 
[0393] Each coefficient can thus be written: 



<* base (S) 
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X* = Xf + aWj + n i5 with 1 < i < N. 



[0394] The detection of the inserted watermarking signal w consists of constructing a hypothesis test for replying to 
the question "has the signal w been inserted in X*?". 

[0395] The advantage of constructing a hypothesis test is that its result makes it possible to give a theoretical confidence 
level to the decision since a hypothesis test is standardised with respect to known distributions. 

[0396] For a class of statistical tests based on correlation, it can be shown that this detection amplitude can be 
calculated from a calculation function of the type: 



noise and a, b and c are constants which depend on the distribution of the watermarking signal w and on the required 
detection probability level p. 

[0397] Thus this minimum value cc D is determined according notably to the length P of the watermarking signal w, a 
detection probability level p, the variance of the coefficients to be modulated, the distribution of the watermarking signal 
and possibly the variance of an additive noise decorrelated from the coefficients modelling any distortions applied to the 
coefficients. 

[0398] Conversely, it is also possible to express the minimum detection length, denoted P°, as a function of the 
weighting coefficient a, for a fixed detection probability level: 



[0399] It would be possible to use this result to check the parameters which make it possible to insert a watermarking 
signal w which will be robust to any additive noise with an energy less than or equal to a n 2 . 

[0400] As illustrated on the curves in Figure 14, from the visual amplitude oc v calculated as a function of the length P 
of the watermarking signal to be inserted, the insertion method includes, in general terms, astep of determining compatible 
values of the length P of the watermarking signal w, of the weighting coefficient a and of the probability level p so that 
the value of the weighting coefficient determined a is less than or equal to the maximum value or visual amplitude oc v 
and greater than or equal to the minimum value or detection amplitude oc D for these determined values of the length P 
of the watermarking signal w and the detection probability level p. 

[0401 ] The weighting coefficient a, the length P of the watermarking signal w and the detection probability level p must 
thus be determined so that the following inequality is satisfied: 



[0402] As illustrated in Figure 14, in order to calculate the minimum value P° of the watermarking signal w for which 
the aforementioned inequality will be satisfied, the formula linking the length P of the watermarking signal w to the 
amplitude of the weighting coefficient a is used, with a = a^N) and with a detection probability level of p%. 
[0403] In this example a standardised statistical test is taken, based on the correlation between the watermarked 
signal X* and the inserted watermarking signal w, as proposed for example in the document entitled "On Resolving 
Rightful Ownerships of Digital Images by Invisible Watermarks" by W. ZENG and B. LIU, in Proc. ICI P 97, pages 552-555, 
October 1997. 

[0404] It can be shown that, for this test: 




yjbP-C 



2 2 

where <7 x is equal to the variance of the coefficients to be modulated, <J n is equal to the variance of the additive 




CCD < Ct <OCv 
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[0405] A calculation step E604 thus makes it possible to calculate the three aforementioned constants a, b and c so that 



b = E 2 (w 2 ), 

and 



[0406] The values E(w 2 ) and cr 2 2 can be calculated theoretically as a function of the distribution of w. For the uniform 
distribution considered here, the theoretical values are E(i^)= 1/3 and (j\ = 4/45. 

[0407] The value t 0 is a threshold value, dependent on the probability of correct detection p : t 0 = 2 x t p . The values 
of tp are percentile of a centred Gaussian law and have a variance equal to 1 . These values are tabulated in the state 
of the art, for example in the book entitled "Probability and Statistics" by A. PAPOULIS, Prentice-Hall, 1990. 
[0408] For example, tp = 3,291 corresponds to p = 0,9995, that is to say to a detection probability level equal to 99,95%. 
[0409] In this numerical example it is considered that any additive noise is not taken into account so that a n = 0. 
[0410] In an estimation step E605, the variance of the set of coefficients X is also estimated. The variance is estimated 
according to the following formula: 



<rl=—1b( x i- M xY 

with 



i .v 

M x = — Tx t 
x n£ 

[0411] Then, for a detection probability fixed at 99.95%, the function linking the length P of the watermarking signal 
w to be inserted and the weighting coefficient a is then derived therefrom: 

p _ 14,44 x<t 2 + 3,84a 2 
0,11a 2 

[0412] At the calculation step E603, the minimum length P° is determined from the above formula for a weighting 
coefficient equal to the maximum value oc v (N) calculated at the calculation step E602. 

[0413] The insertion method includes, in this example embodiment, a step E606 of comparing this calculated length 
P° with the cardinal number N of the set of modulatable coefficients X. 
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[0414] If, at this comparison step E606, the calculated length P° is less than or equal to the cardinal number N of the 
set of modulatable coefficients X, it is deduced therefrom that the insertion of a watermarking signal w of length P° is 
possible and meets the required imperceptibility and detection criteria. 

[0415] In reality, it is possible to effect an insertion of a watermarking signal w which meets these criteria since the 
5 length P of this watermarking signal lies between the calculated value P° and the cardinal number N of the set of 
modulatable coefficients X. 

[0416] Then the insertion proper of a watermarking signal of length P is carried out by modulating a subset of coefficients 
with a cardinal number P according to the linear modulation level, with a weighting coefficient a equal to the calculated 
maximum value oc v (N). 

10 [0417] The insertion method includes, in the usual fashion, a reverse spatio-frequency transformation step which 
makes it possible to obtain, at the output of the coder, the watermarked image I*. 

[041 8] Where insertion is possible, it is also possible to calculate the number of theoretical bits which can be inserted. 
This is because, at the decoder, the decoding operation consists of detecting whether a pseudo-random sequence w 
has been inserted. The response being binary (yes/no), it can be considered that a pseudo-random sequence of length 
15 p makes it possible to insert and extract an information bit in the image. By repeating the insertion operation on several 
subsets of spectral coefficients, it is thus possible to insert and extract several information bits. 

[041 9] Consequently, since in this example there are a set of coefficients of cardinal number N and since only a number 
P° of coefficients are necessary to insert and detect an information bit, it is possible to derive therefrom the total number 
Q of information bits which can be inserted, referred to as the capacity of the image: Q = N/P°. 
20 [0420] Moreover, if at the comparison step E606 the calculated length P° is greater than the cardinal number N of the 
set of modulatable coefficients X, it is a priori deduced therefrom that the insertion of a watermarking signal meeting the 
criteria of imperceptibility and detectability is not possible. 

[0421] In this case, the insertion method can if necessary include an additional step of calculating the detection 
probability level p' obtained for a minimum value oc D of the weighting coefficient equal to the maximum value 0Cy(N) 
25 calculated for a length P of the watermarking signal w equal to the cardinal number N of the set of modulatable coefficients 
X. 

[0422] If this detection probability level, which will be less than that fixed initially at 99.95%, is judged satisfactory, the 
insertion of a watermarking signal of length equal to N can be effected, that is to say all the coefficients of the sub-band 
of frequency HH 1 are modulated. 
30 [0423] In a second embodiment of the invention, the insertion method can include a step of calculating a threshold 
value P* of the length of the watermarking signal w, as illustrated in Figure 14. 

[0424] This threshold value P* is determined so that the minimum value oc D of the weighting coefficient is equal to the 
maximum value oc v of the weighting coefficient for a length P of the watermarking signal w equal to this threshold value 
P* and a predetermined detection probability level p. 
35 [0425] This calculation step makes it possible to determine, from formulae expressing the detection amplitude oc D and 
the visual amplitude oc v notably as a function of the length P of the watermarking signal w, the minimum value P* of the 
insertion signal from which the weighting coefficient a can satisfy the inequality oc D < a <oc v . 

[0426] In Figure 14, this minimum value P* is given by the intersection of the two curves representing the variation in 
the visual amplitude and the variation in the detection amplitude as a function of the length P of the inserted water- 
^o marking signal w, for a predetermined detection probability level. 

[0427] Thus, in general terms, the insertion method according to the invention has the advantage of being able to 
control a priori and theoretically the detection probability and imperceptibility of the watermarking signal which it is sought 
to insert in an image. 

[0428] Naturally, many modifications could be made to the example embodiments described above without departing 
45 from the scope of the invention. 

[0429] Thus the insertion technique used could be applied to the raw digital image, without undergoing any spatio- 
frequency transformation prior to the modulation of the coefficients. 

[0430] I n this case, the modulated coefficients are coefficients representing the digital image solely in the spatial domain. 
[0431] Moreover, the spatio-frequency transformation applied to the image can use analysis and synthesis filters other 
50 than those described above, or even be a transformation other than discrete Fourier transformation, discrete cosine 
transformation or Fourier-Mellin transformation. These transformations are used currently in conventional processings 
of digital images. 

[0432] In addition, the invention is not limited to the example embodiments described but also concerns any embodiment 
where the values of the weighting coefficient a, of the length P of the watermarking signal and of the detection probability 
55 level p are determined in a correlated fashion so as to satisfy the inequality oc D < a < a v . 

[0433] The present invention also concerns a device for determining a partitioning and a decoding device which are 
adapted to implement the partitioning and decoding methods described previously. 

[0434] These devices for determining a partitioning and decoding can be used in a computer 50 as illustrated in Figure 
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15, independently or in the same computer 50. 

[0435] The partitioning determination device comprises partitioning means 500, 502, 503 adapted to perform an 
adaptive partitioning according to a criterion of detectability of an information bit inserted on each region. 
[0436] These partitioning means use on coding a maximum value oc v of the weighting coefficient ensuring impercep- 
5 tibility of the watermarking signal for modulating all the coefficients representing said digital image, possibly with a 
modulation factor k, dependent on each coefficient, as described previously. 

[0437] This device for determining a partitioning also has means 500, 502, 503 of comparing the size P of each region 
of the partitioning with a minimum size P stat corresponding to minimum size of a statistically significant sample for the 
detection of a watermarking signal inserted in said region. 

10 [0438] The device also has means 500, 502, 503 for applying apre-distortion before calculating a detectability criterion. 
[0439] The partitioning means comprise means 500, 502, 503 of modulating all the coefficients of each region by 
inserting the watermarking signal, means 500, 502, 503 of calculating a detectability amplitude from the detectability 
criterion and means 500, 502, 503 of validating the detectability criterion, adapted to compare the detectability amplitude 
with a predetermined threshold value Tc. 

15 [0440] The partitioning means 500, 502, 503 are adapted to effect a partitioning by iterative division of the numerical 
image and, for each region at a rank in the said partitioning, the validation means 500, 502, 503 are adapted to validate 
a partitioning of a rank higher in the region if and only if there exist at least two subregions of the region for which the 
detectability criterion is validated. 

[0441 ] This device also comprises means 500, 502, 503 of incrementing the value of a capacity Q of the digital image, 
20 adapted, for each region at a rank in the partitioning, to increment the capacity value Q when the partitioning of rank 
directly higher is not validated, the capacity Q being equal to the cardinal number of the distinct regions obtained by the 
partitioning on which the detectability criterion is validated. 

[0442] In another embodiment, this device for determining a partitioning comprises incrementation means 500, 502, 
503 adapted, for each region R at a rank in the partitioning, to increment the value of a capacity Q when the partitioning 
25 of rank directly higher is not validated, comparison means 500, 502, 503 adapted to compare the value of the capacity 
Q with a fixed capacity value Qs and reiteration means 500, 502, 503 adapted to reiterate the partitioning for a region 
having the highest detectability amongst all the other regions to be processed when the capacity value Q is less than 
the fixed capacity value Qs. 

[0443] In a third embodiment of the invention, the partitioning means 500, 502, 503 are adapted to effect a partitioning 
30 by iterative merging of the digital image and there are adapted, for each region of the partitioning for which the detectability 
criterion is not validated, to merge the region with another region of the partitioning, and preferably with another region 
of the partitioning for which the detectability criterion is not validated. 

[0444] In an entirely similar fashion, the device for decoding a watermarking signal in a digital image comprises means 
500, 502, 503 of partitioning the digital image to be decoded into distinct regions, adapted to perform an adaptive 
35 partitioning as a function of a criterion of detectability of an information bit obtained by demodulating the coefficients on 
each region. 

[0445] It also has means 500, 502, 503 of comparing the size P of each region R of the partitioning with a minimum 
size P stat corresponding to the minimum size of a statistically significant sample for the detection of a watermarking 
signal inserted in the region. 

40 [0446] It comprises means 500, 502, 503 of calculating a detection amplitude from the detectability criterion and means 
500, 502, 503 of validating the detectability criterion, adapted to compare the detection amplitude with a predetermined 
threshold value Td for decoding. 

[0447] This predetermined threshold value Td for decoding is less than the predetermined threshold value Tc used 
during the partitioning determination method described previously. 
45 [0448] In a preferred embodiment, the device for determining a partition includes means 500, 502, 503 for checking 
a criterion of detectability, over each region, of an information bit inserted by means of a pseudo-random sequence with 
average of zero by modulation of the coefficients of the region. 

[0449] It also includes means 500, 502, 503 for generating a centred pseudo-random sequence equal in size to the 
cardinal number of the set of coefficients representative of an image, here equal to N. 
50 [0450] Means 500, 502, 503 of modulating the set of coefficients X by the centred pseudo-random sequence in order 
to insert the same information bit onto this set of coefficients are also incorporated into the device for determining a 
partition 1 1 as shown on figure 1 . 

[0451] The latter may also includes means 500, 502, 503 for applying a distortion D which makes it possible to provide 
for and guarantee the robustness of the inserted watermark in the case of certain types of distortion, typically a com- 
55 pression of the digital image in order to store it. 

[0452] These means of applying a distortion make it possible actually to simulate a predetermined distortion over the 
set of modulated coefficients. Here they comprise, whenever the coefficients of the set X are spectral coefficients, means 
500, 502, 503 for inverse spectral transformation of the set of coefficients modulated in order to obtain a watermarked 
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image, means 500, 502, 503 for applying a distortion to this watermarked image, and means 500, 502, 503 for spectral 
transformation of the watermarked image in order to recover the set of modulated coefficients after distortion. 
[0453] In this example, the inverse-spectral-transformation means are suitable for implementing multi-resolution spec- 
tral wavelet recomposition and the forward spectral transformation means are suitable for again reiterating the breakdown 
of the digital image I into wavelets. 

[0454] The device for determining a partition 11 , in a first way, also includes means for partitioning by iterative com- 
bination, on the basis of an initial segmentation of the set of coefficients X. 

[0455] For preference, this initial segmentation makes it possible to segment the set of coefficients X into a set of 
regions of a size which is statistically significant for validating the detectability criterion over these regions. 
[0456] The means of partitioning 500, 502, 503 by iterative combination are suitable, for each region of the partition 
for which the detectability criterion is not satisfied, for combining the region with another region of the partition. 
[0457] Obviously, the partitioning means 500, 502, 503 could also be suitable for forming this partition of the image 
by iterative division of the digital image according to asecond way of performing the preferred embodiment of the invention. 
[0458] This partition-determining device 1 1 makes it possible, as output, to obtain a partition into regions of the set of 
coefficients X, with it being possible to insert an information bit onto each of these regions with a reliable detection 
probability and robustness to a given distortion. The maximum number of regions which can be obtained by this adaptive 
partitioning of the set X corresponds to the practical insertion capacity of this set X. 

[0459] As illustrated in Figure 1 , an insertion device 1 0, incorporating the device for determining a partition 1 1 described 
above, moreover comprises means 13 for reading a watermark to be inserted, containing several information bits. By 
way of an example which is not at all limiting, a binary watermark may be 1 0001 1 1 , representing, for example, the coding 
of the name of the author of the image to be watermarked. 

[0460] The insertion means 1 2 proper may be conventional and make it possible, by modulation, to insert the various 
information bits of the watermarkonto the various insertion carriers determined by the device for determining apartition 1 1 . 
[0461] According to the invention, the device for inserting a watermarking signal w comprises means 500, 502, 503 
of spatio-frequency transformation of an image I, and for example analysis filters associated with decimators by two 
adapted to effect a wavelet decomposition of an image I. It also has reverse spatio-frequency recomposition means 500, 
502, 503 for recomposing the image I after the insertion of the watermarking signal in the domain transformed into sub- 
bands. 

[0462] It also comprises means 500, 502, 503 of inserting a watermarking signal w of length P adapted to modulate 
a subset of coefficients of cardinal number P according to the linear model described previously using a weighting 
coefficient a. 

[0463] This insertion device also comprises means 500, 502, 503 of calculating the maximum value a v of the weighting 
coefficient as a function of the length P of the watermarking signal w ensuring imperceptibility of the watermarking signal. 
These calculation means use the function oc^P, T) described above and dependent on the length P of the watermarking 
signal w and the spectral transformation T used for decomposing the image signal I. 

[0464] In general terms, this insertion device also has means 500, 502, 503 of determining compatible values of the 
length P of the watermarking signal w, the weighting coefficient a and the probability level p so that the value of the 
weighting coefficient determined is less than or equal to the maximum value oc v and greater than or equal to the minimum 
value cc D for the determined values of the length P of the watermarking signal w and the detection probability level p. 
[0465] The calculation means 500, 502, 503 are also adapted to determine a minimum value or detection amplitude 
oc D using the calculation functioning described above of the type: 



[0466] In one embodiment of the invention, the insertion device also has means of calculating 500, 502, 503 the length 
P° of the watermarking signal w for a minimum value oc D of the weighting coefficient equal to a maximum value oc^N) 
calculated for a length P of the watermarking signals w equal to the cardinal number N of the set of modulatable coefficients 
ensuring imperceptibility of the watermarking signal w at the limit of perceptibility and for a predetermined detection 
probability level p; and means 500, 502, 503 of comparing this length P° calculated with the cardinal number N of the 
set of modulatable coefficients X. 

[0467] It also comprises means 500, 502, 503 of calculating a detection probability level p' obtained for a minimum 
value oc D of the weighting coefficient equal to the maximum value Oy calculated for a length P of the watermarking signal 
equal to the cardinal number N of the set of modulatable coefficients X. 

[0468] In another embodiment, the insertion device comprises means 500, 502, 503 of calculating a threshold value 
P* of the length of the watermarking signal w determined so that the minimum value oc D of the weighting coefficient is 



cc D {P) = 




■JbP-c 
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equal to the maximum value cc v of the weighting coefficient of the length P of the watermarking signal w equal to this 
threshold value P* and a predetermined detection probability level p. 

[0469] All the previously described devices are incorporated in a microprocessor 500 of the computer 50, a read-only 
memory 502 being adapted to store a program for determining a partitioning and/or inserting a watermarking signal 
and/or decoding the watermarking signal and a random access memory 503 comprising registers adapted to store 
variables modified during the running of the programs. 

[0470] The microprocessor 500 is integrated into the computer 50, which can be connected to different peripherals, 
for example a digital camera 507 or a microphone 51 1 , by means of an input/output card 506 in order to receive and 
store documents. 

[0471] The digital camera 507 makes it possible notably to supply images to be authenticated by inserting a water- 
marking signal. 

[0472] This computer 50 has a communication interface 512 connected to a communication network 513 in order to 
receive, if necessary, images to be watermarked. 

[0473] The computer 50 also has document storage means, such as a hard disk 508, or is adapted to cooperate, by 
means of a disk drive 509, with removable document storage means such as diskettes 51 0. 

[0474] These fixed or removable storage means can also contain the code of the different methods according to the 
invention which, once read by the microprocessor 500, will be stored on the hard disk 508. 

[0475] By way of variant, the program enabling the devices to implement the invention can be stored in the read only 
memory 502. 

[0476] As a second variant, the program can be received in order to be stored as described above by means of the 
communication network 51 3. 

[0477] The computer 50 also has a screen 504 for serving, for example, as an interface with an operator by means 
of the keyboard 514 or any other means. 

[0478] The central unit 500 (CPU) will execute the instructions relating to the implementation of the invention. On 
powering up, the programs and methods relating to the invention stored with a non-volatile memory, for example the 
read only memory 502, are transferred into the random access memory 503 (RAM) which will then contain the executable 
code of the invention and the variables necessary for implementing the invention. 

[0479] This random access memory 503 contains different registers for storing the variables necessary for the running 
of the program, and notably a register for storing the coefficients of the regions R which are temporarily modulated or 
demodulated in order to verify the detectability criterion, the minimum value P stat , and the threshold value Tc and Td 
respectively on coding and decoding. 

[0480] This random-access memory 503 also includes an X h Xj register for storing the coefficients of the regions which 
are modulated or demodulated temporarily so as to check the detectability criterion and a register w for storing the 
centred pseudo-random sequence generated in order to determine the partition of the set X. 

[0481] Finally, the random access memory 503 contains a register for storing the table of visual amplitudes of cc base 
(T), a register for storing the spectral coefficients X, a register for calculating the visual amplitude oc^N) at the limit of 
imperceptibility, a register for storing the values of the constants a, b, c and of the variance o x 2 necessary for determining 
the detection amplitude cc D as a function of the length P of the watermarking signal w and a register for storing the 
minimum value P° calculated. 

[0482] A communication bus 501 affords communication between the different sub-elements of the computer 50 or 
linked to it. The representation of the bus 501 is not limitative and notably the microprocessor 500 is able to indicate 
instructions to any sub-element directly or by means of another sub-element. 



Claims 

1. Method of inserting a watermarking signal W into a digital image I by modulating coefficients X of distinct regions 
R representing said digital image I, comprising a step of determining a partitioning into distinct regions of the digital 
image characterised in that the partitioning into regions is effected by an adaptive partitioning as a function of a 
criterion of statistical detectability of an information bit inserted on each region R. 

2. Method according to Claim 1 , wherein the modulating of coefficients X is function of a weighting coefficient, char- 
acterised in that the detectability criterion is calculated by using a maximum value ocv of the weighting coefficient 
ensuring imperceptibility of the watermarking signal for modulating all the coefficients representing said digital image. 

3. Method of determining a partitioning according to claim 2, characterised in that a maximum value (ocv) can be 
calculated as a function of the length P of the watermarking signal w by: 
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a y (P,S,w) = 



where 0Cbase( s ) is a base value, dependent on the transformation used S and on the (base) sub-band considered 
for the insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of said 
sub-band is modulated, 
P is strictly greater than 2, and 

E[|iv|P] is the mathematical expectation of the function \w\P. 

4. Method according to one of Claims 2 or 3, characterised in that the detectability criterion is calculated by using, 
for each coefficient (Xj) to be modulated, a weighting coefficient in accordance with a law of the form otj = kj.oc v , 
where kj is a modulation factor dependent on the coefficients situated close to the relevant coefficient (Xj) on said 
region and Oy is equal to said maximum value of the weighting coefficient. 

5. Method according to one of claims 1 to 4, characterised in that it comprises a step (E307) of applying a distortion 
before calculating the detectability criterion. 

6. Method according to one of Claims 1 to 5, characterised in that it includes a prior step (E105) of comparing the 
size P of each region R of the partitioning with a minimum size P stat corresponding to the minimum size with a 
statistically significant sample for the detection of a watermarking signal inserted in said region. 

7. Method according to one of Claims 1 to 6, characterised in that, for each region R of the partitioning, all the 
coefficients of said region are modulated by inserting the watermarking signal W, a detectability amplitude T is 
calculated from said detectability criterion and said detectability criterion is validated by comparing said detectability 
amplitude T with a predetermined threshold value Tc. 

8. Method according to claim 7, characterised in that said detectability criterion is validated if a minimum modulation 
length P min (x) calculated for ensuring detectability of the inserted watermarking signal W is less than or equal to the 
cardinal number of said region R. 

9. Method according to claim 1 , characterised in that it comprises the following steps: 

- generation (E501 ) of a centred pseudo-random sequence w equal in size to the cardinal number N of a set of 
coefficients representative of at least a part of an image I, formed from centred pseudo-random sub-sequences 
w k ; 

- modulation (E502) of the said set of coefficients by the said centred pseudo-random sequence w in order to 
insert the same information bit on the said set of coefficients; and 

- checking (E506, E507, E514, E516), over each region, of a criterion of detectability of the information bit 
inserted by modulation. 

10. Method according to Claim 9, characterised in that it comprises a single step of applying a distortion (E503) to 
thesaidsetof modulated coefficients, before the step of checking (E506, E507, E514, E51 6) of a detectability criterion. 

1 1 . Method according to Claim 1 0, characterised in that the said set of coefficients is a set of spectral coefficients X 
originating from a spectral transformation of a digital image, and in that the said step (E503) of applying a distortion 
includes the following sub-steps: 

- inverse spectral transformation (E503a) of the said set of modulated coefficients in order to obtain a water- 
marked image; 

- application (E503b) of a distortion to the said watermarked image; and 

- spectral transformation (E503c) of the watermarked image in order to obtain a set of modulated coefficients 
after distortion. 

12. Method according to one of Claims 9 to 11 , characterised in that the centred pseudo-random sequence w is formed 
from pseudo-random sub-sequences w k of identical size. 
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13. Method according to one of Claims 9 to 12, characterised in that the pseudo-random sub-sequences w k have a 
size greater than or equal to a minimum number of coefficients forming a region of the said partition of a size which 
is statistically significant for validating the said criterion of detectability over the said region. 

14. Method according to one of Claims 9 to 13, characterised in that, at the generation step (E501), each centred 
pseudo-random sub-sequence w k is created by generating half of the pseudo-random samples of the said sub- 
sequence via a pseudo-random number generator and by symmetrising the said pseudo-random samples generated. 

15. Method according to Claim 14, characterised in that, at the generation step (E501), each centred pseudo-random 
sub-sequence w k is created by moreover permutating the said pseudo-random samples using a secret key. 

16. Method according to one of Claims 9 to 13, characterised in that, at the generation step (E501), each centred 
pseudo-random sub-sequence w k is created by generating pseudo-random random samples via a pseudo-random 
number generator and by redistributing over each sample the sum of the said pseudo-random samples generated. 

17. Method according to one of Claims 7 to 16, characterised in that a partitioning by iterative division (E103) of the 
digital image is effected and in that, for each region at a rank of said partitioning, a partitioning of a directly higher 
rank of said region is validated if and only if there exist at least two sub-regions of said region for which said 
detectability criterion is validated. 

18. Method according to Claim 17, characterised in that, for each region at a rank of said partitioning, the value of a 
capacity Q of the digital image is incremented (E1 12) when the partitioning of directly higher rank is not validated, 
said capacity Q being equal to the cardinal number of the distinct regions obtained by said partitioning on which the 
detectability criterion is validated. 

19. Method according to Claim 18, characterised in that, for each region at a rank of said partitioning, the value of a 
capacity Q is incremented (E1 12) when the partitioning of immediately high rank is not validated, said value of the 
capacity Q is compared (E1 15) with a fixed capacity value Qs and the partitioning is reiterated for a region having 
the highest detectability amplitude amongst all the other regions to be processed when said capacity value is lower 
than said fixed capacity value. 

20. Method according to one of Claims 7 to 16, characterised in that a partitioning by iterative merging of the digital 
image is effected and in that, for each region of said partitioning for which said detectability criterion is not validated, 
said region is merged with another region of said partitioning. 

21. Method according to Claim 20, characterised in that said other region of the partitioning is if possible a region for 
which said detectability criterion is not validated. 

22. Method of decoding a watermarking signal in a digital image, inserted by modulating coefficients X representing 
said image in distinct regions forming a partitioning of said digital image, characterised in that a partitioning into 
distinct regions R of the digital image to be decoded is determined by effecting an adaptive partitioning as a function 
of a criterion of statistical detectability of an information bit obtained by demodulating coefficients on each region R. 

23. Decoding method according to Claim 22, characterised in that it includes a prior step (E201) of comparing the 
size P of each region of the partitioning with a minimum size P stat corresponding to the minimum size of a statistically 
significant sample for the detection of a watermarking signal inserted in said region. 

24. Decoding method according to one of Claims 22 or 23, characterised in that, for each region R of the partitioning, 
a detection amplitude T is calculated (E203) from said detectability criterion and said detectability criterion is validated 
(E203) by comparing said detection amplitude T with a predetermined threshold value Td for decoding. 

25. Decoding method according to Claim 24, characterised in that the predetermined threshold value Td for decoding 
is less than said predetermined threshold value Tc used during the inserting method according to Claim 7. 

26. Device for inserting a watermarking signal into a digital image by modulating coefficients of regions representing 
said digital image, comprising means for partitioning into distinct regions the digital image characterised in that 
partitioning means (500, 502, 503) are adapted to effect an adaptive partitioning as a function of a criterion of 
statistical detectability of an information bit inserted on each region. 
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27. Device according to claim 26, characterised in that it also includes means (500, 502, 503) for applying a distortion 
before calculating the detectability criterion. 

28. Device according to one of claims 26 to 27, characterised in that it also has means (500, 502, 503) of comparing 
the size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistically 
significant sample for the detection of a watermarking signal inserted in the said region. 

29. Device according to one of Claims 26 to 28, characterised in that the partitioning means (500, 502, 503) comprise 
means (500, 502, 503) of modulating all the coefficients of each region by inserting the watermarking signal, means 
(500, 502, 503) of calculating a detectability amplitude from said detectability criterion and means (500, 502, 503) 
of validating said detectability criterion, adapted to compare said detectability amplitude with a predetermined thresh- 
old value. 

30. Device according to claim 26, characterised in that it includes: 

- means (500, 502, 503) for generating a centred pseudo-random sequence w equal in size to the cardinal 
number N of a set of coefficients representative of at least a part of an image I, formed by centred pseudo- 
random sub-sequences w k ; 

- means (500, 502, 503) for modulation of the said set of coefficients by the said centred pseudo-random 
sequence w in order to insert the same information bit on the said set of coefficients; and 

- means (500, 502, 503) for checking, over each region, a criterion of detectability of the information bit inserted 
by modulation. 

31. Device according to Claim 30, characterised in that it comprises means (500, 502, 503) for applying a distortion 
to the said set of modulated coefficients. 

32. Device according to Claim 31 , characterised in that the said set of coefficients is a set of spectral coefficients Xi 
originating from a spectral transformation of a digital image, and in that the said means for applying a distortion 
comprise: 

- means (500, 502, 503) for inverse spectral transformation of the said set of modulated coefficients in order to 
obtain a watermarked image; 

- means (500, 502, 503) for application of a distortion to the said watermarked image; and 

- means (500, 502, 503) for spectral transformation of the watermarked image in order to obtain aset of modulated 
coefficients after distortion. 

33. Device according to one of Claims 29 to 32, characterised in that the partitioning means (500, 502, 503) are 
adapted to effect a partitioning by iterative division of the digital image and in that, for each region at a rank of said 
partitioning, the validation means are adapted to validate a partitioning of a directly higher rank of said region if and 
only if there exist at least two sub-regions of said regions for which said detectability criterion is validated. 

34. Device according to Claim 33, characterised in that it also comprises means (500, 502, 503) of incrementing the 
value of a capacity of the digital image adapted, for each region at a rank of said partitioning, to increment the said 
capacity value when the partitioning of directly higher rank is not validated, said capacity being equal to the cardinal 
number of the distinct regions obtained by said partitioning on which the detectability criterion is validated. 

35. Device according to Claim 33, characterised in that it also comprises incrementation means (500, 502, 503) 
adapted, for each region at a rank of said partitioning, to increment the value of a capacity when the partitioning of 
directly higher rank is not validated, comparison means (500, 502, 503) adapted to compare said value of the 
capacity with a fixed capacity value and reiteration means (500, 502, 503) adapted to reiterate the partitioning for 
a region having the highest detectability amplitude amongst all the other regions to be processed when said capacity 
value is lower than said fixed capacity value. 

36. Device according to one of Claims 29 to 32, characterised in that the partitioning means (500, 502, 503) are 
adapted to effect a partitioning by iterative merging of the digital image and are adapted, for each region of said 
partitioning for which said detectability criterion is not validated, to merge said region with another region of said 
partitioning. 
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37. Device according to Claim 36, characterised in that said other region of the partitioning is if possible a region for 
which said detectability criterion is not validated. 

38. Device according to one of Claims 26 to 37, characterised in that it is incorporated in a microprocessor (500), a 
read only memory (502) being adapted to store a program for determining a partitioning and a random access 
memory (503) comprising registers adapted to store variables modified during the running of said program. 

39. Device for decoding a watermarking signal in a digital image, inserted by modulation of coefficients representing 
said image in distinct regions forming a partitioning of said digital image, characterised in that it comprises means 
(500, 502, 503) of partitioning of the digital image to be decoded into distinct regions, adapted to effect an adaptive 
partitioning as a function of a criterion of statistical detectability of an information bit obtained by demodulation of 
the coefficients on each region. 

40. Decoding device according to Claim 39, characterised in that it also has means (500, 502, 503) of comparing the 
size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistical significant 
sample for the detection of a watermarking signal inserted in said region. 

41. Decoding device according to one of Claims 39 or 40, characterised in that it comprises means (500, 502, 503) 
of calculating a detection amplitude from said detectability criterion and means (500, 502, 503) of validating said 
detectability criterion adapted to compare said detection amplitude with a threshold value predetermined for the 
decoding. 

42. Decoding device according to Claim 41 , characterised in that the predetermined threshold value for the decoding 
is less than said predetermined threshold value used during the inserting method according to Claim 7. 

43. Decoding device according to one of Claims 39 to 42, characterised in that it is incorporated in a microprocessor 
(500), a readonly memory (502) being adapted to store a program for decoding a watermarking signal and a random 
access memory (503) comprising registers adapted to store variables modified during the running of said program. 

44. Computer, characterised in that it comprises means adapted to implement an inserting method according to one 
of Claims 1 to 21 . 

45. Computer, characterised in that it comprises means adapted to implement the decoding method according to one 
of Claims 22 to 25. 

46. Information storage means readable by a microprocessor, characterised in that it stores a program implementing 
an inserting method according to one of Claims 1 to 21 . 

47. Information storage means readable by a microprocessor, characterised in that it stores a program implementing 
a decoding method according to one of Claims 22 to 25. 



Patentanspruche 

1. Verfahren zum Einfugen eines Wasserzeichensignals W in ein digitales Bild I durch Modulieren von Koeffizienten 
X von verschiedenen Bereichen R, die das digitale Bild I darstellen, mit einem Schritt zum Bestimmen einer Parti- 
tionierung des digitalen Bildes in verschiedene Bereiche, dadurch gekennzeichnet, dass die Partitionierung in 
Bereiche durch eine adaptive Partitionierung als Funktion eines Kriteriums einer statistischen Erfassbarkeit eines 
in jedem Bereich R eingefugten Informationsbits ausgefuhrt wird. 

2. Verfahren gemaB Anspruch 1, bei dem das Modulieren von Koeffizienten X eine Funktion eines Gewichtungskoef- 
fizienten ist, dadurch gekennzeichnet, dass das Erfassbarkeitskriterium durch Verwendung eines Hochstwerts 
ocv des Gewichtungskoeffizienten berechnet wird, der eine Nichtwahrnehmbarkeit des Wasserzeichensignals fur 
eine Modulation aller Koeffizienten sicherstellt, die das digitale Bild darstellen. 

3. Verfahren zum Bestimmen einer Partitionierung gemaG Anspruch 2, dadurch gekennzeichnet, dass ein Hochst- 
wert (ocv) als Funktion der Lange P des Wasserzeichensignals w berechnet werden kann durch: 
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a v (P r S,W)= (X Basis( S ^ , 



wobei oc Basis (S) ein Basiswert, der von der verwendeten Transformation S und von dem fur die Einfugung in Betracht 
gezogenen (Basis-)Teilband abhangig ist, von dem maximalen Gewichtungskoeffizienten ist, der eine Nichtwahr- 
10 nehmbarkeit sicherstellt, wenn ein einziger Koeffizient des Teilbands moduliert wird, 

P groBer ist als 2, und 

E[|w|P] der mathematische Erwartungswert der Funktion |w|P ist. 

4. Verfahren gemaB einem der Anspruche 2 oder 3, dadurch gekennzeichnet, dass das Erfassbarkeitskriterium 
15 berechnet wird, indem fur jeden zu modulierenden Koeffizienten (Xj) ein Gewichtungskoeffizient gemaB einem 

Gesetz in der Form otj = k h a v verwendet wird, wobei kj ein Modulationsfaktor ist, der von den Koeffizienten abhangig 
ist, die sich nahe an dem maBgeblichen Koeffizienten (Xj) in dem Bereich befinden, und Oy gleich dem Hochstwert 
des Gewichtungskoeffizienten ist. 

20 5. Verfahren gemaB einem der Anspruche 1 bis 4, dadurch gekennzeichnet, dass es einen Schritt (E307) zum 
Anwenden einer Verzerrung vor einer Berechnung des Erfassbarkeitskriteriums aufweist. 

6. Verfahren gemaB einem der Anspruche 1 bis 5, dadurch gekennzeichnet, dass es einen vorherigen Schritt (E1 05) 
zum Vergleichen der GroBe P von jedem Bereich R der Partitionierung mit einer MindestgroBe P stat aufweist, die 

25 der MindestgroBe eines statistisch signifikanten Abtastwerts fur die Erfassung eines in dem Bereich eingefugten 

Wasserzeichensignals entspricht. 

7. Verfahren gemaB einem der Anspruche 1 bis 6, dadurch gekennzeichnet, dass fur jeden Bereich R der Partitio- 
nierung alle Koeffizienten des Bereichs moduliert werden, indem das Wasserzeichensignal W eingefugt wird, eine 

30 Erfassbarkeitsamplitude T aus dem Erfassbarkeitskriterium berechnet wird und das Erfassbarkeitskriterium durch 

Vergleich der Erfassbarkeitsamplitude T mit einem vorbestimmten Schwellenwert Tc validiert wird. 

8. Verfahren gemaB Anspruch 7, dadurch gekennzeichnet, dass das Erfassbarkeitskriterium validiert wird, wenn 
eine minimale Modulationslange P m j n (x), die zum Sicherstellen einer Erfassbarkeitdes eingefugten Wasserzeichen- 

35 signals W berechnet wird, kleiner oder gleich der Kardinalzahl des Bereichs R ist. 

9. Verfahren gemaB Anspruch 1 , dadurch gekennzeichnet, dass es die folgenden Schritte aufweist: 

- Erzeugung (E501) einer zentrierten Pseudozufallssequenz w, die groBenmaBig gleich der Kardinalzahl N 
^o einer Menge von Koeffizienten ist, die fur zumindest einen Teil eines Bildes I reprasentativ sind, wobei sie aus 

zentrierten Pseudozufallsuntersequenzen w k gebildet wird; 

- Modulation (E502) der Menge von Koeffizienten durch die zentrierte Pseudozufallssequenz w, urn das gleiche 
Informationsbit in die Menge von Koeffizienten einzufugen; und 

- Uberprufung (E506, E507, E514, E516) eines Kriteriums einer Erfassbarkeitdes durch Modulation eingefugten 
45 Informationsbits uber jeden Bereich. 

10. Verfahren gemaB Anspruch 9, dadurch gekennzeichnet, dass es einen einzigen Schritt zum Anwenden einer 
Verzerrung (E503) auf die Menge von modulierten Koeffizienten vor dem Schritt der Uberprufung (E506, E507, 
E514, E516) eines Erfassbarkeitskriteriums aufweist. 

50 

11. Verfahren gemaB Anspruch 10, dadurch gekennzeichnet, dass die Menge von Koeffizienten eine Menge von 
Spektralkoeffizienten X ist, die aus einer Spektraltransformation eines digitalen Bildes stammen, und dass der Schritt 
(E503) zum Anwenden einer Verzerrung die folgenden Unterschritte aufweist: 

55 - inverse Spektraltransformation (E503a) der Menge von modulierten Koeffizienten, urn ein mit Wasserzeichen 

versehenes Bild zu erhalten; 

- Anwendung (E503b) einer Verzerrung auf das mit Wasserzeichen versehene Bild; und 

- Spektraltransformation (E503c) des mit Wasserzeichen versehenen Bildes, urn eine Menge von modulierten 
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Koeffizienten nach einer Verzerrung zu erhalten. 

12. Verfahren gemaG einem der Anspruche 9 bis 11 , dadurch gekennzeichnet, dass die zentrierte Pseudozufallsse- 
quenz w aus Pseudozufallsuntersequenzen w k identischer GroGe gebildet wird. 

13. Verfahren gemaG einem der Anspruche 9 bis 12, dadurch gekennzeichnet, dass die Pseudozufallsuntersequen- 
zen w k eine GroGe aufweisen, die groGer oder gleich einer minimalen Anzahl von Koeffizienten ist, die einen Bereich 
der Partitionierung mit einer GroGe bilden, die fur eine Validierung des Erfassbarkeitskriteriums uber den Bereich 
statistisch signifikant ist. 

14. Verfahren gemaG einem der Anspruche 9 bis 13, dadurch gekennzeichnet, dass im Erzeugungsschritt (E501) 
jede zentrierte Pseudozufallsuntersequenz w k erstellt wird, indem die Halfte der Pseudozufallsabtastwerte der Un- 
tersequenz uber einen Pseudozufallszahlengenerator erzeugt wird, und indem die erzeugten Pseudozufallsabtast- 
werte symmetrisch gemacht werden. 

15. Verfahren gemaG Anspruch 1 4, dadurch gekennzeichnet, dass im Erzeugungsschritt (E501 ) jede zentrierte Pseu- 
dozufallsuntersequenz w k erstellt wird, indem auGerdem die Pseudozufallsabtastwerte unter Verwendung eines 
geheimen Schlussels permutiert werden. 

16. Verfahren gemaG einem der Anspruche 9 bis 13, dadurch gekennzeichnet, dass im Erzeugungsschritt (E501) 
jede zentrierte Pseudozufallsuntersequenz w k erstellt wird, indem Pseudozufallsabtastwerte uber einen Pseudozu- 
fallszahlengenerator erzeugt werden, und indem die Summe der erzeugten Pseudozufallsabtastwerte uber jeden 
Abtastwert umverteilt wird. 

17. Verfahren gemaG einem der Anspruche 7 bis 1 6, dadurch gekennzeichnet, dass eine Partitionierung durch iterative 
Aufteilung (E1 03)desdigitalen Bildes ausgefuhrt wird, und dass fur jeden Bereich auf einer Ebeneder Partitionierung 
eine Partitionierung einer direkt hoheren Ebene des Bereichs validiert wird, wenn und nur wenn zumindest zwei 
Unterbereiche des Bereichs existieren, fur die das Erfassbarkeitskriterium validiert wird. 

18. Verfahren gemaG Anspruch 17, dadurch gekennzeichnet, dass fur jeden Bereich auf einer Ebene der Partitio- 
nierung der Wert einer Kapazitat Q des digitalen Bildes inkrementiert wird (E1 12), wenn die Partitionierung einer 
direkt hoheren Ebene nicht validiert wird, wobei die Kapazitat Q gleich der Kardinalzahl der verschiedenen durch 
die Partitionierung erhaltenen Bereiche ist, in denen das Erfassbarkeitskriterium validiert wird. 

19. Verfahren gemaG Anspruch 18, dadurch gekennzeichnet, dass fur jeden Bereich auf einer Ebene der Partitio- 
nierung der Wert einer Kapazitat Q inkrementiert wird (E1 12), wenn die Partitionierung einer unmittelbar hoheren 
Ebene nicht validiert wird, wobei der Wert der Kapazitat Q mit einem festen Kapazitatswert Qs verglichen wird 
(E115) und die Partitionierung fur einen Bereich wiederholt wird, der unter all den anderen zu verarbeitenden 
Bereichen die hochste Erfassbarkeitsamplitude aufweist, wenn der Kapazitatswert kleiner ist als der feste Kapazi- 
tatswert. 

20. Verfahren gemaG einem der Anspruche 7 bis 16, dadurch gekennzeichnet, dass eine Partitionierung durch ite- 
ratives Vereinigen des digitalen Bildes ausgefuhrt wird, und dass fur jeden Bereich der Partitionierung, fur den das 
Erfassbarkeitskriterium nicht validiert wird, dieser Bereich mit einem anderen Bereich der Partitionierung vereinigt 
wird. 

21. Verfahren gemaG Anspruch 20, dadurch gekennzeichnet, dass der andere Bereich der Partitionierung, wenn 
moglich, ein Bereich ist, fur den das Erfassbarkeitskriterium nicht validiert wird. 

22. Verfahren zum Decodieren eines Wasserzeichensignals in einem digitalen Bild, das durch Modulieren von Koeffi- 
zienten X, die das Bild darstellen, in verschiedenen Bereichen, die eine Partitionierung des digitalen Bildes bilden, 
eingefugt ist, dadurch gekennzeichnet, dass eine Partitionierung in verschiedene Bereiche R des zu decodieren- 
den digitalen Bildes bestimmt wird, indem eine adaptive Partitionierung als Funktion eines Kriteriums einer statisti- 
schen Erfassbarkeit eines Informationsbits ausgefuhrt wird, das durch Demodulieren von Koeffizienten in jedem 
Bereich R erhalten wird. 

23. Decodierverfahren gemaG Anspruch 22, dadurch gekennzeichnet, dass es einen vorherigen Schritt (E201 ) zum 
Vergleichen der GroGe P von jedem Bereich der Partitionierung mit einer MindestgroGe P stat aufweist, die der 
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MindestgroGe eines statistisch signifikanten Abtastwerts fur die Erfassung eines in dem Bereich eingefugten Was- 
serzeichensignals entspricht. 

24. Decodierverfahren gemaG einem der Anspruche 22 Oder 23, dadurch gekennzeichnet, dass fur jeden Bereich R 
der Partitionierung eine Erfassungsamplitude T aus dem Erfassbarkeitskriterium berechnet wird (E203) und das 
Erfassbarkeitskriterium validiert wird (E203), indem die Erfassungsamplitude T mit einem vorbestimmten Schwel- 
lenwert Td fur eine Decodierung verglichen wird. 

25. Decodierverfahren gemaG Anspruch 24, dadurch gekennzeichnet, dass der vorbestimmten Schwellenwert Td 
fur eine Decodierung kleiner ist als der vorbestimmte Schwellenwert Tc, der wahrend des Einfugungsverfahrens 
gemaG Anspruch 7 verwendet wird. 

26. Vorrichtung zum Einfugen eines Wasserzeichensignals in ein digitales Bild durch Modulieren von Koeffizienten von 
Bereichen, die das digitale Bild darstellen, mit einer Einrichtung zum Partitionieren des digitalen Bildes in verschie- 
dene Bereiche, dadurch gekennzeichnet, dass eine Partitioniereinrichtung (500, 502, 503) angepasst ist zum 
Ausfuhren einer adaptiven Partitionierung als Funktion eines Kriteriums einer statistischen Erfassbarkeit eines in 
jedem Bereich eingefugten Informationsbits. 

27. Vorrichtung gemaG Anspruch 26, dadurch gekennzeichnet, dass sie auch eine Einrichtung (500, 502, 503) zum 
Anwenden einer Verzerrung vor einer Berechnung des Erfassbarkeitskriteriums umfasst. 

28. Vorrichtung gemaG einem der Anspruche 26 bis 27, dadurch gekennzeichnet, dass sie auch eine Einrichtung 
(500, 502, 503) zum Vergleichen der GroGe von jedem Bereich der Partitionierung mit einer MindestgroGe umfasst, 
die der MindestgroGe eines statistisch signifikanten Abtastwerts fur die Erfassung eines in dem Bereich eingefugten 
Wasserzeichensignals entspricht. 

29. Vorrichtung gemaG einem der Anspruche 26 bis 28, dadurch gekennzeichnet, dass die Partitioniereinrichtung 
(500, 502, 503) eine Einrichtung (500, 502, 503) zum Modulieren aller Koeffizienten von jedem Bereich durch 
Einfugen des Wasserzeichensignals, eine Einrichtung (500, 502, 503) zum Berechnen einer Erfassbarkeitsamplitude 
aus dem Erfassbarkeitskriterium und eine Einrichtung (500, 502, 503) zum Validieren des Erfassbarkeitskriteriums 
aufweist, die angepasst ist zum Vergleichen der Erfassbarkeitsamplitude mit einem vorbestimmten Schwellenwert. 

30. Vorrichtung gemaG Anspruch 26, dadurch gekennzeichnet, dass sie umfasst: 

- eine Einrichtung (500, 502, 503) zum Erzeugen einer zentrierten Pseudozufallssequenz w, die groGenmaGig 
gleich der Kardinalzahl N einer Menge von Koeffizienten ist, die fur zumindest einen Teil eines Bildes I repra- 
sentativ sind, wobei sie durch zentrierte Pseudozufallsuntersequenzen w k gebildet wird; 

- eine Einrichtung (500, 502, 503) fur eine Modulation der Menge von Koeffizienten durch die zentrierte Pseu- 
dozufallssequenz w, urn das gleiche Informationsbit in die Menge von Koeffizienten einzufugen; und 

- eine Einrichtung (500, 502, 503) zum Uberprufen eines Kriteriums einer Erfassbarkeit des durch Modulation 
eingefugten Informationsbits Ciber jeden Bereich. 

31. Vorrichtung gemaG Anspruch 30, dadurch gekennzeichnet, dass sie eine Einrichtung (500, 502, 503) zum An- 
wenden einer Verzerrung auf die Menge von modulierten Koeffizienten aufweist. 

32. Vorrichtung gemaG Anspruch 31, dadurch gekennzeichnet, dass die Menge von Koeffizienten ein Menge von 
Spektralkoeffizienten Xi ist, die aus einer Spektraltransformation eines digitalen Bildes stammen, und dass die 
Einrichtung zum Anwenden einer Verzerrung aufweist: 

- eine Einrichtung (500, 502, 503) fur eine inverse Spektraltransformation der Menge von modulierten Koeffi- 
zienten, urn ein mit Wasserzeichen versehenes Bild zu erhalten; 

- eine Einrichtung (500, 502, 503) zur Anwendung einer Verzerrung auf das mit Wasserzeichen versehene Bild; 
und 

- eine Einrichtung (500, 502, 503) fur eine Spektraltransformation des mit Wasserzeichen versehenen Bildes, 
urn eine Menge von modulierten Koeffizienten nach einer Verzerrung zu erhalten. 

33. Vorrichtung gemaG einem der Anspruche 29 bis 32, dadurch gekennzeichnet, dass die Partitioniereinrichtung 
(500, 502, 503) angepasst ist zum Ausfuhren einer Partitionierung durch iterative Aufteilung des digitalen Bildes, 
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und dass die Validierungseinrichtung fur jeden Bereich auf einer Ebene der Partitionierung angepasst ist, eine 
Partitionierung einer direkt hoheren Ebene des Bereichs zu validieren, wenn und nur wenn zumindest zwei Unter- 
bereiche des Bereichs existieren, fur die das Erfassbarkeitskriterium validiert wird. 

34. Vorrichtung gemaG Anspruch 33, dadurch gekennzeichnet, dass sie auch eine Einrichtung (500, 502, 503) zum 
Inkrementieren des Werts einer Kapazitat des digitalen Bildes aufweist, die fur jeden Bereich auf einer Ebene der 
Partitionierung angepasst ist zum Inkrementieren des Kapazitats werts, wenn die Partitionierung einer direkt hoheren 
Ebene nicht validiert wird, wobei die Kapazitat gleich der Kardinalzahl der verschiedenen durch die Partitionierung 
erhaltenen Bereiche ist, in denen das Erfassbarkeitskriterium validiert wird. 

35. Vorrichtung gemaG Anspruch 33, dadurch gekennzeichnet, dass sie auch aufweist: eine Inkrementierungsein- 
richtung (500, 502, 503), die fur jeden Bereich auf einer Ebene der Partitionierung angepasst ist zum Inkrementieren 
des Werts einer Kapazitat, wenn die Partitionierung einer direkt hoheren Ebene nicht validiert wird, eine Vergleichs- 
einrichtung (500, 502, 503), die angepasst ist zum Vergleichen des Werts der Kapazitat mit einem festen Kapazi- 
tatswert, und eine Wiederholungseinrichtung (500, 502, 503), die angepasst ist zum Wiederholen der Partitionierung 
fur einen Bereich, der unter all den anderen zu verarbeitenden Bereichen die hochste Erfassbarkeitsamplitude 
aufweist, wenn der Kapazitatswert kleiner ist als der feste Kapazitatswert. 

36. Vorrichtung gemaG einem der Anspruche 29 bis 32, dadurch gekennzeichnet, dass die Partitioniereinrichtung 
(500, 502, 503) angepasst ist zum Ausfuhren einer Partitionierung durch iteratives Vereinigen des digitalen Bildes 
und fur jeden Bereich der Partitionierung, fur den das Erfassbarkeitskriterium nicht validiert wird, angepasst ist zum 
Vereinigen des Bereichs mit einem anderen Bereich der Partitionierung. 

37. Vorrichtung gemaG Anspruch 36, dadurch gekennzeichnet, dass der andere Bereich der Partitionierung, wenn 
mdglich, ein Bereich ist, fur den das Erfassbarkeitskriterium nicht validiert wird. 

38. Vorrichtung gemaG einem der Anspruche 26 bis 37, dadurch gekennzeichnet, dass sie in einem Mikroprozessor 
(500) eingebunden ist, wobei ein Festwertspeicher (502) angepasst ist zum Speichern eines Programms zum 
Bestimmen einer Partitionierung und ein Direktzugriffsspeicher (503) Register aufweist, die angepasst sind zum 
Speichern von Variablen, die wahrend des Ablaufs des Programms modifiziert werden. 

39. Vorrichtung zum Decodieren eines Wasserzeichensignals in einem digitalen Bild, das durch Modulieren von Koef- 
fizienten, die das Bild darstellen, in verschiedenen Bereichen, die eine Partitionierung des digitalen Bildes bilden, 
eingefugt ist, dadurch gekennzeichnet, dass sie eine Einrichtung (500, 502, 503) zum Partitionieren des zu 
decodierenden digitalen Bildes in verschiedene Bereiche aufweist, die angepasst ist zum Ausfuhren einer adaptiven 
Partitionierung als Funktion eines Kriteriums einer statistischen Erfassbarkeit eines Informationsbits, das durch 
Demodulieren der Koeffizienten in jedem Bereich erhalten wird. 

40. Decodiervorrichtung gemaG Anspruch 39, dadurch gekennzeichnet, dass sie auch eine Einrichtung (500, 502, 
503) zum Vergleichen der GroGe von jedem Bereich der Partitionierung mit einer MindestgroGe aufweist, die der 
MindestgroGe eines statistisch signifikanten Abtastwerts fur die Erfassung eines in dem Bereich eingefugten Was- 
serzeichensignals entspricht. 

41 . Decodiervorrichtung gemaG einem der Anspruche 39 oder 40, dadurch gekennzeichnet, dass sie eine Einrichtung 
(500, 502, 503) zum Berechnen einer Erfassungsamplitude aus dem Erfassbarkeitskriterium und eine Einrichtung 
(500, 502, 503) zum Validieren des Erfassbarkeitskriteriums aufweist, die angepasst ist zum Vergleichen der Er- 
fassungsamplitude mit einem fur die Decodierung vorbestimmten Schwellenwert. 

42. Decodiervorrichtung gemaG Anspruch 41, dadurch gekennzeichnet, dass der vorbestimmte Schwellenwert fur 
die Decodierung kleiner ist als der vorbestimmte Schwellenwert, der wahrend des Einfugungsverfahrens gemaG 
Anspruch 7 verwendet wird. 

43. Decodiervorrichtung gemaG einem der Anspruche 39 bis 42, dadurch gekennzeichnet, dass sie in einem Mikro- 
prozessor (500) eingebunden ist, wobei ein Festwertspeicher (502) angepasst ist zum Speichern eines Programms 
zum Decodieren eines Wasserzeichensignals und ein Direktzugriffsspeicher (503) Register aufweist, die angepasst 
sind zum Speichern von Variablen, die wahrend des Ablaufs des Programms modifiziert werden. 

44. Computer, dadurch gekennzeichnet, dass dieser eine Einrichtung aufweist, die angepasst ist zum Implementieren 
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eines Einfugungsverfahrens gemaB einem der Anspruche 1 bis 21 . 

45. Computer, dadurch gekennzeichnet, dass dieser eine Einrichtung aufweist, die angepasst ist zum Implementieren 
des Decodierverfahrens gemaB einem der Anspruche 22 bis 25. 

46. Informationsspeichereinrichtung, die von einem Mikroprozessor lesbar ist, dadurch gekennzeichnet, dass sie ein 
Programm speichert, das ein Einfugungsverfahren gemaB einem der Anspruche 1 bis 21 implementiert. 

47. Informationsspeichereinrichtung, die von einem Mikroprozessor lesbar ist, dadurch gekennzeichnet, dass sie ein 
Programm speichert, das ein Decodierverfahren gemaB einem der Anspruche 22 bis 25 implementiert. 



Revendications 

1. Procede d'insertion d'un signal de marquage W dans une image numerique I en modulant des coefficients X de 
regions distinctes R representant ladite image numerique I, comportant une etape consistant a determiner un par- 
titionnement en regions distinctes de I'image numerique caracterise en ce que le partitionnement en regions est 
effectue par I'intermediaire d'un partitionnement adaptatif en fonction d'un critere de detectabilite statistique d'un 
bit d'informations insere sur chaque region R. 

2. Procede selon la revendication 1 , dans lequel la modulation de coefficients X estfonction d'un coefficient ponderateur, 
caracterise en ce que le critere de detectabilite est calcule en utilisant une valeur maximum ocv du coefficient 
ponderateur assurant I'imperceptibilite du signal de marquage pour moduler tous les coefficients representant ladite 
image numerique. 

3. Procede de determination d'un partitionnement selon la revendication 2, caracterise en ce qu'une valeur maximum 
(ocv) peut etre calculee en fonction de la longueur P du signal de marquage w par : 



ou oCbase(S) estune valeur de base, dependantede la transformation utilisee S etde lasous-bande (base) consideree 
pour I'insertion, du coefficient ponderateur maximum assurant I'imperceptibilite lorsqu'un coefficient unique de ladite 
sous-bande est module, 
P est strictement superieur a 2, et 

E[|w|P] est I'esperance mathematique de la fonction |w|P. 

4. Procede selon I'une des revendications 2 et 3, caracterise en ce que le critere de detectabilite est calcule en 
utilisant, pour chaque coefficient (X|) a moduler, un coefficient ponderateur conformement a une loi de la forme oc, 
= kj.a v , ou kj est un facteur de modulation dependant des coefficients situes a proximite du coefficient approprie 
(Xj) sur ladite region et est egal a ladite valeur maximum du coefficient ponderateur. 

5. Procede selon I'une des revendications 1 a 4, caracterise en ce qu'il comporte une etape (E307) consistant a 
appliquer une distorsion avant de calculer le critere de detectabilite. 

6. Procede selon I'une des revendications 1 a 5, caracterise en ce qu'il inclut une etape prealable (E1 05) consistant 
a comparer la taille P de chaque region R du partitionnement a une taille minimum P stat correspondant a la taille 
minimum avec un echantillon statistiquement significatif pour la detection d'un signal de marquage insere dans 
ladite region. 

7. Procede selon I'une des revendications 1 a 6, caracterise en ce que, pour chaque region R du partitionnement, 
tous les coefficients de ladite region sont modules en inserant le signal de marquage W, une amplitude de detectabilite 
T est calculee a partir dudit critere de detectabilite et ledit critere de detectabilite est valide en comparant ladite 
amplitude de detectabilite T a une valeur de seuil predeterminee Tc. 
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8. Procede selon la revendication 7, caracterise en ce que ledit critere de detectabilite est valide si une longueur de 
modulation minimum P m j n (X) calculee pour assurer la detectabilite du signal de marquage insere W est inferieure 
ou egale au nombre cardinal de ladite region R. 

5 9. Procede selon la revendication 1 , caracterise en ce qu'il comporte les etapes suivantes : 

- la generation (E501) d'une sequence pseudo-aleatoire centree w egale en taille au nombre cardinal N d'un 
ensemble de coefficients representatifs d'au moins une partie d'une image I, formee a partir de sous-sequences 
pseudo-aleatoires centrees w k , 

10 - la modulation (E502) dudit ensemble de coefficients par ladite sequence pseudo-aleatoire centree w afin 

d'inserer le meme bit d'informations sur ledit ensemble de coefficients, et 

- le controle (E506, E507, E514, E516), sur chaque region, d'un critere de detectabilite du bit d'informations 
insere par modulation. 

15 1 o. Procede selon la revendication 9, caracterise en ce qu'il comporte une etape unique d'application d'une distorsion 
(E503) audit ensemble de coefficients modules, avant I'etape de controle (E506, E507, E514, E516) d'un critere de 
detectabilite. 

11. Procede selon la revendication 10, caracterise en ce que ledit ensemble de coefficients est un ensemble de 
20 coefficients spectraux X provenant d'une transformation spectrale d'une image numerique, et en ce que ladite etape 

(E503) consistant a appliquer une distorsion inclut les sous-etapes suivantes : 

- la transformation spectrale inverse (E503a) dudit ensemble de coefficients modules afin d'obtenir une image 
marquee, 

25 - 1'application (E503b) d'une distorsion a ladite image marquee, et 

- la transformation spectrale (E503c) de I'image marquee afin d'obtenir un ensemble de coefficients modules 
apres distorsion. 

12. Procede selon Tune des revendications 9 a 11 , caracterise en ce que la sequence pseudo-aleatoire centree w est 
30 formee a partir de sous-sequences pseudo-aleatoires w k de taille identique. 

13. Procede selon I'une des revendications 9 a 12, caracterise en ce que les sous-sequences pseudo-aleatoires w k 
ont une taille superieure ou egale a un nombre minimum de coefficients formant une region de ladite partition d'une 
taille qui est statistiquement significative pour valider ledit critere de detectabilite sur ladite region. 

35 

14. Procede selon I'une des revendications 9 a 13, caracterise en ce que, a I'etape de generation (E501), chaque 
sous-sequence pseudo-aleatoire centree w k est creee en generant la moitie des echantillons pseudo-aleatoires de 
ladite sous-sequence via un generateur de nombres pseudo-aleatoires et en symetrisant lesdits echantillons pseudo- 
aleatoires generes. 

40 

15. Procede selon la revendication 1 4, caracterise en ce que, a I'etape de generation (E501 ), chaque sous-sequence 
pseudo-aleatoire centree w k est creee en permutant en outre lesdits echantillons pseudo-aleatoires en utilisant une 
cle secrete. 

^5 16. Procede selon I'une des revendications 9 a 13, caracterise en ce que, a I'etape de generation (E501), chaque 
sous-sequence pseudo-aleatoire centree w k est creee en generant des echantillons aleatoires pseudo-aleatoires 
viaun generateur de nombres pseudo-aleatoires eten redistribuant sur chaque echantillon lasomme desdits echan- 
tillons pseudo-aleatoires generes. 

50 17. Procede selon I'une des revendications. 7 a 16, caracterise en ce qu'un partitionnement par division iterative 
(E1 03) de I'image numerique esteffectue eten ceque, pour chaque region au niveau d'un rang dudit partitionnement, 
un partitionnement d'un rang directement superieur de ladite region est valide si et seulement s'il existe au moins 
deux sous-regions de ladite region pour lesquelles ledit critere de detectabilite est valide. 

55 18. Procede selon la revendication 17, caracterise en ce que, pour chaque region au niveau d'un rang dudit partition- 
nement, la valeur d'une capacite Q de I'image numerique est incrementee (E112) lorsque le partitionnement de 
rang directement superieur n'est pas valide, ladite capacite Q etant egale au nombre cardinal des regions distinctes 
obtenues par ledit partitionnement sur lesquelles le critere de detectabilite est valide. 
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19. Procede selon la revendication 18, caracterise en ce que, pour chaque region au niveau d'un rang dudit partition- 
nement, la valeur d'une capacite Q est incrementee (E112) lorsque le partitionnement du rang immediatement 
superieur n'est pas valide, ladite valeur de la capacite Q est comparee (E1 15) a une valeur de capacite fixe Qs et 
le partitionnement est repete pour une region ayant ramplitude de detectabilite la plus elevee parmi toutes les autres 
regions a traiter lorsque ladite valeur de capacite est inferieure a ladite valeur de capacite fixe. 

20. Procede selon Tune des revendications 7 a 16, caracterise en ce qu'un partitionnement par fusion iterative de 
I'image numerique est effectue et en ce que, pour chaque region dudit partitionnement pour laquelle ledit critere 
de detectabilite n'est pas valide, ladite region est fusionnee avec une autre region dudit partitionnement. 

21. Procede selon la revendication 20, caracterise en ce que ladite autre region du partitionnement est si possible 
une region pour laquelle ledit critere de detectabilite n'est pas valide. 

22. Procede de decodage d'un signal de marquage dans une image numerique, insere en modulant des coefficients X 
representant ladite image dans des regions distinctes formant un partitionnement de ladite image numerique, ca- 
racterise en ce qu'un partitionnement en regions distinctes R de I'image numerique a decoder est determine en 
effectuant un partitionnement adaptatif en fonction d'un critere de detectabilite statistique d'un bit d'informations 
obtenu en demodulant des coefficients sur chaque region R. 

23. Procede de decodage selon la revendication 22, caracterise en ce qu'il inclut une etape prealable (E201 ) consistant 
a comparer la taille P de chaque region du partitionnement a une taille minimum P stat correspondant a la taille 
minimum d'un echantillon statistiquement significatif pour la detection d'un signal de marquage insere dans ladite 
region. 

24. Procede de decodage conformement a Tune des revendications 22 et 23, caracterise en ce que, pour chaque 
region R du partitionnement, une amplitude de detection T est calculee (E203) a partir dudit critere de detectabilite 
et ledit critere de detectabilite est valide (E203) en comparant ladite amplitude de detection T a une valeur de seuil 
predeterminee Td en vue d'un decodage. 

25. Procede de decodage selon la revendication 24, caracterise en ce que la valeur de seuil predeterminee Td pour 
effectuer un decodage est inferieure a ladite valeur de seuil predeterminee Tc utilisee pendant le procede d'insertion 
conformement a la revendication 7. 

26. Dispositif d'insertion d'un signal de marquage dans une image numerique en modulant des coefficients de regions 
representant ladite image numerique, comportant des moyens pour partitionner en regions distinctes I'image nu- 
merique caracterise en ce que des moyens de partitionnement (500, 502, 503) sont adaptes pour effectuer un 
partitionnement adaptatif en fonction d'un critere de detectabilite statistique d'un bit d'informations insere sur chaque 
region. 

27. Dispositif selon la revendication 26, caracterise en ce qu'il inclut egalement des moyens (500, 502, 503) pour 
appliquer une distorsion avant de calculer le critere de detectabilite. 

28. Dispositif selon I'une des revendications 26 et 27, caracterise en ce qu'il dispose egalement de moyens (500, 502, 
503) de comparaison de la taille de chaque region du partitionnement a une taille minimum correspondant a la taille 
minimum d'un echantillon statistiquement significatif pour la detection d'un signal de marquage insere dans ladite 
region. 

29. Dispositif selon I'une des revendications 26 a 28, caracterise en ce que les moyens de partitionnement (500, 502, 
503) comportent des moyens (500, 502, 503) de modulation de tous les coefficients de chaque region en inserant 
le signal de marquage, des moyens (500, 502, 503) de calcul d'une amplitude de detectabilite a partir dudit critere 
de detectabilite et des moyens (500, 502, 503) de validation dudit critere de detectabilite, adaptes pour comparer 
ladite amplitude de detectabilite a une valeur de seuil predeterminee. 

30. Dispositif selon la revendication 26, caracterise en ce qu'il inclut : 

- des moyens (500, 502, 503) pour generer une sequence pseudo-aleatoire centree w egale en taille au nombre 
cardinal N d'un ensemble de coefficients representatifs d'au moins une partie d'une image I, formee par des 
sous-sequences pseudo-aleatoires centrees w k , 
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- des moyens (500, 502, 503) pour moduler ledit ensemble de coefficients par ladite sequence pseudo-aleatoire 
centree w afin d'inserer le meme bit d'informations sur ledit ensemble de coefficients, et 

- des moyens (500, 502, 503) pour controler, sur chaque region, un critere de detectabilite du bit d'informations 
insere par modulation. 

31. Dispositif selon la revendication 30, caracterise en ce qu'il comporte des moyens (500, 502, 503) pour appliquer 
une distorsion audit ensemble de coefficients modules. 

32. Dispositif selon la revendication 31, caracterise en ce que ledit ensemble de coefficients est un ensemble de 
coefficients spectraux Xi provenant d'une transformation spectrale d'une image numerique, et en ce que lesdits 
moyens destines a appliquer une distorsion comportent : 

- des moyens (500, 502, 503) pour une transformation spectrale inverse dudit ensemble de coefficients modules 
afin d'obtenir une image marquee, 

- des moyens (500, 502, 503) pour I'application d'une distorsion a ladite image marquee, et 

- des moyens (500, 502, 503) pour une transformation spectrale de I'image marquee afin d'obtenir un ensemble 
de coefficients modules apres la distorsion. 

33. Dispositif selon I'une des revendications 29 a 32, caracterise en ce que les moyens de partitionnement (500, 502, 
503) sont adaptes pour effectuer un partitionnement par division iterative de I'image numerique et en ce que, pour 
chaque region au niveau d'un rang dudit partitionnement, les moyens de validation sont adaptes pour valider un 
partitionnement d'un rang directement superieur de ladite region si et seulement s'il existe au moins deux sous- 
regions desdites regions pour lesquelles ledit critere de detectabilite est valide. 

34. Dispositif selon la revendication 33, caracterise en ce qu'il comporte egalement des moyens (500, 502, 503) 
decrementation de la valeur d'une capacite de I'image numerique adaptee, pour chaque region au niveau d'un 
rang dudit partitionnement, pour incrementer ladite valeur de capacite lorsque le partitionnement du rang directement 
superieur n'est pas valide, ladite capacite etant egale au nombre cardinal des regions distinctes obtenues par ledit 
partitionnement sur lequel le critere de detectabilite est valide. 

35. Dispositif selon la revendication 33, caracterise en ce qu'il comporte egalement des moyens decrementation 
(500, 502, 503) adaptes, pour chaque region au niveau d'un rang dudit partitionnement, pour incrementer la valeur 
d'une capacite lorsque le partitionnement de rang directement superieur n'est pas valide, des moyens de compa- 
raison (500, 502, 503) adaptes pour comparer ladite valeur de la capacite a une valeur de capacite fixe et des 
moyens de repetition (500, 502, 503) adaptes pour repeter le partitionnement pour une region ayant I'amplitude de 
detectabilite la plus elevee parmi toutes les autres regions a traiter lorsque ladite valeur de capacite est inferieure 
a ladite valeur de capacite fixe. 

36. Dispositif selon I'une des revendications 29 a 32, caracterise en ce que les moyens de partitionnement (500, 502, 
503) sont adaptes pour effectuer un partitionnement par fusion iterative de I'image numerique et sont adaptes, pour 
chaque region dudit partitionnement pour laquelle ledit critere de detectabilite n'est pas valide, pour fusionner ladite 
region avec une autre region dudit partitionnement. 

37. Dispositif selon la revendication 36, caracterise en ce que ladite autre region du partitionnement est si possible 
une region pour laquelle ledit critere de detectabilite n'est pas valide. 

38. Dispositif selon I'une des revendications 26 a 37, caracterise en ce qu'il est incorpore dans un microprocesseur 
(500), une memoire a lecture seule (502) etant adaptee pour memoriser un programme afin de determiner un 
partitionnement et une memoire a acces direct (503) comportant des registres adaptes pour memoriser des variables 
modifiees pendant I'execution dudit programme. 

39. Dispositif de decodage d'un signal de marquage dans une image numerique, insere par modulation de coefficients 
representant ladite image dans des regions distinctes formant un partitionnement de ladite image numerique, ca- 
racterise en ce qu'il comporte des moyens (500, 502, 503) de partitionnement de I'image numerique a decoder 
en regions distinctes, adaptes pour effectuer un partitionnement adaptatif en fonction d'un critere de detectabilite 
statistique d'un bit d'informations obtenu par demodulation des coefficients sur chaque region. 

40. Dispositif de decodage selon la revendication 39, caracterise en ce qu'il dispose egalement de moyens (500, 502, 
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503) de comparaison de la taille de chaque region du partitionnement a une taille minimum correspondant a la taille 
minimum d'un echantillon statistique significatif pour la detection d'un signal de marquage insere dans ladite region. 

41 . Dispositif de decodage selon Tune des revendications 39 et 40, caracterise en ce qu'il comporte des moyens (500, 
502, 503) de calcul d'une amplitude de detection a partir dudit critere de detectabilite et des moyens (500, 502, 503) 
de validation dudit critere de detectabilite adaptes pour comparer ladite amplitude de detection a une valeur de seuil 
predeterminee pour le decodage. 

42. Dispositif de decodage selon la revendication 41 , caracterise en ce que la valeur de seuil predeterminee pour le 
decodage est inferieure a ladite valeur de seuil predeterminee utilisee pendant le procede d'insertion selon la 
revendication 7. 

43. Dispositif de decodage selon Tune des revendications 39 a 42, caracterise en ce qu'il est incorpore dans un 
microprocesseur (500), une memoire a lecture seule (502) etant adaptee pour memoriser un programme pour 
decoder un signal de marquage et une memoire a acces direct (503) comportant des registres adaptes pour me- 
moriser des variables modifiees pendant I'execution dudit programme. 

44. Ordinateur, caracterise en ce qu'il comporte des moyens adaptes pour mettre en oeuvre un procede d'insertion 
conformement a Tune des revendications 1 a 21 . 

45. Ordinateur, caracterise en ce qu'il comporte des moyens adaptes pour mettre en oeuvre le procede de decodage 
selon Tune des revendications 22 a 25. 

46. Moyens de memorisation d'informations lisibles par un microprocesseur, caracterises en ce qu'ils memorisent un 
programme mettant en oeuvre un procede d'insertion conformement a I'une des revendications 1 a 21 . 

47. Moyens de memorisation d'informations lisibles par un microprocesseur, caracterises en ce qu'ils memorisent un 
programme mettant en oeuvre un procede de decodage conformement a I'une des revendications 22 a 25. 
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